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lar dimensional reduction, in the context of finite temperature field theory. We first 
briefly review the formalism of thermal field theory and show how dimensional re- 
duction emerges as the high-temperature limit for static quantities. Then we apply 
dimensional reduction to two distinct problems, the pressure of electroweak theory 
and the screening masses of mesonic operators in hot QCD, and point out the simi- 
larities. We summarize the results and discuss their validity, while leaving all details 
to original research articles. 
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Chapter 1 

Introduction 



The physics of interactions between elementary particles is described to an amazing ac- 
curacy by the standard model of particle physics. It ties three of the four fundamental 
interactions, namely the electromagnetic, weak and strong interactions, together under 
the conceptual framework of relativistic quantum field theory. Scattering processes and 
bound states involving few particles are well described by the model, although many open 
questions, mostly related to strongly interacting states, remain. In the energy region cur- 
rently accessible to experiments we have therefore full reason to believe that this theory 
is correct. 

When matter is heated high above everyday temperatures, its neutral constituents are 
torn apart into an interacting plasma of elementary particles. At temperatures of the same 
order or higher than the particle masses this necessitates combining quantum statistical 
mechanics with relativistic field theory. The interactions between individual particles are 
still governed by the standard model interactions, but the effects of hot medium change 
their long-distance behavior and give rise to many-particle collective modes. 

Experimentally such extreme conditions are accessible in relativistic heavy ion collisions 
currently produced at the Relativistic Heavy Ion Collider (RHIC) in Brookhaven, and, 
starting this year, also at Large Hadron Collider (LHC) at CERN. In these experiments two 
heavy nuclei collide against each other, forming a finite volume of extremely hot matter. 
The matter described by the theory of strong interactions, quantum chromodynamics 
(QCD), goes through a phase transition to a deconfined phase of color-charged particles 
forming a quark-gluon plasma, which then rapidly cools as it expands. This kind of 
temperatures were also present in the very early universe, whose expansion is sensitive to 
the equation of state of both QCD and electroweak matter. 

The formalism for finite temperature quantum field theory arises naturally from the path 
integral quantization of field theories. The time coordinate is extended to complex values 
to account for varying the fields over statistical ensemble, and the functional integral is over 
all field configurations periodic or antiperiodic in the imaginary time. When temperature 
is larger than any other scale in the process, the excitations in the imaginary time can 
be integrated out and the physics of static quantities is described by a three-dimensional 
effective theory. This is known as dimensional reduction [5]. The effective theory can be 
systematically derived, and it exhibits the same infrared behavior as the full theory. At 
finite temperatures the main advantage in using a dimensionally reduced effective theory 
in perturbative computations is the ability to systematically treat the various infrared 
divergences, as well as the resummations needed to cure them, in a simpler setting. 
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Dimensional reduction has been successfully applied over the years to compute many 
bosonic quantities both perturbatively and in combination with lattice simulations. In the 
QCD sector, the three-dimensional formulation known as EQCD has made it possible to 
perturbatively compute the pressure up to the last perturbative order g^ln g [6-13], and 
the result has also been extended to nonzero chemical potentials [14]. Lattice implementa- 
tions of EQCD have been used to compute the static correlation lengths of various gluonic 
operators [15-20]. There are also recent developments in formulating an effective theory 
preserving the spontaneously broken Z{3) symmetry of the deconfined phase [21], which 
is explicitly broken in EQCD [22]. Besides QCD, the electroweak symmetry breaking has 
also been solved in detail using lattice simulations in dimensionally reduced effective the- 
ory [23-28] , motivated by the possibility of a first order electroweak phase transition being 
the origin of the observed baryon asymmetry in the universe. 

There are only few applications of dimensional reduction to fermionic observables, be- 
cause the fermion fields are integrated out from the three-dimensional effective theory. This 
simplifies the computation of bosonic quantities tremendously, but the accessible fermionic 
observables are then limited to those that can be inferred from vacuum or bosonic ones, 
such as quark number susceptibilities Xij = d'^p/dfiidfij [29]. Systematic application of di- 
mensional reduction to fermionic operators was developed in [30], inspired by the progress 
in heavy quark effective theories. 

The use of dimensional reduction is restricted to time-independent quantities. It should 
be mentioned here that for real-time computations there exists another scheme of resum- 
ming the light particle self-energy corrections to regulate some of the infrared divergences, 
namely the hard thermal loop (HTL) approximation [31,32]. Both schemes succesfully 
resum the one-loop infrared divergences, but in general the HTL Green's functions are 
more complicated, since they carry the full analytic structure of the original theory. It 
is also very hard to systematically improve the HTL approximation beyond the leading 
order. 

In this thesis we study two applications of dimensional reduction to the standard model, 
the perturbative evaluation of the electroweak pressure and the next-to-leading order cor- 
rection to screening masses of mesonic operators. The thesis is organized as follows. In 
chapter [2] we first review the formalism of thermal quantum field theory, and then dis- 
cuss dimensional reduction in the context of general effective theories in section 12.21 In 
chapter [3] we go through the computation of the electroweak theory pressure, with special 
attention given to the behavior near the phase transition. We combine the result with the 
previously known QCD pressure in section [3T3l and study the convergence of the series and 
the deviation from the ideal gas for physical values of parameters. Results for a simpler, 
weakly coupled SU(2) -|- Higgs theory are also shown for comparison. 

In chapter H] we review our work on meson correlators. After a short motivation using 
linear response theory, we compute the leading order correlators at zero and finite density. 
Then we proceed to derive a dimensionally reduced effective theory for the lowest fermionic 
modes and solve the 0{g'^) corrections to screening masses. Finally, we compare with 
recent lattice determinations of the masses and discuss the differences. Chapter [5] contains 
our conclusions. 
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Chapter 2 

Thermal field theory 



In this chapter we will first review how the thermodynamical treatment of quantum field 
theory can be formulated in terms of Euclidean path integrals. We then proceed to discuss 
dimensional reduction, which is the underlying effective theory method used in all the 
research papers included in this thesis. 

2.1 Basic thermodynamics of quantum fields 

The statistical properties of relativistic quantum field theory are most naturally described 
using the grand canonical ensemble. Since particles can be spontaneously created and 
annihilated, the micro canonical or canonical ensembles with fixed particle numbers cannot 
be built, but instead one would have to use the conserved quantities like electric charge. 
To avoid this kind of complicated constraints on field configurations, it is generally easier 
to fix the mean values of energy and conserved commuting number operators using the 
Lagrange multipliers (3 = \/T and /ij, respectively. This is the grand canonical ensemble. 

The thermodynamical properties of the system are given by the partition function and 
its derivatives. In quantum mechanics the partition function is defined as the trace of the 
density matrix p, 

Z(r,y,^i) = Trp = Tre-^^^-'^'^^), (2.1) 

where H and Ni are the Hamiltonian and conserved number operators, respectively. The 
thermal average of an operator is then defined as 

(A) = |TrM, (2.2) 

and the usual thermodynamic quantities like pressure, entropy, energy and particle num- 
bers are given by the partial derivatives 

dT\nZ 
dT 

-pV + TS + fiiNi. (2.3) 

In quantum mechanics the evaluation of the trace in Eq. ()2.ip is simple, one just takes 
any complete orthonormal basis preferably eigenstates olH—fiN if these are known, 

and sums over {n\p\n). The same procedure can in principle be applied to field theory. 



Ni = T— , E 
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where the sum over basis vectors is replaced by a functional integral in the space of field 
configurations. 

Field theories are usually defined in the Lagrangian formalism, and finding the Hamil- 
tonian function required for computation of the partition function in Eq. (12. ip can be 
quite involved, in particular in the context of gauge theories. One has to fix the gauge 
and then carefully separate the canonical variables from auxiliary ones depending on the 
chosen gauge [33]. In addition to the usual canonical equations of motion, the fields are 
constrained by the gauge condition and the field equation for the auxiliary field, which 
can be interpreted as the Gauss' law. 

Once the Hamiltonian has been found, we can insert a complete set of eigenstates 
\<f){x);t) of the field operator (/!>(x) in the Heisenberg picture to compute the partition 
function. This gives 

ZiT,V,f,,) = I m (<A(x);t|e-^(^-^»^')|</>(x);t), (2.4) 

where the integration is over all canonical variables. From the time-dependence of the 
field operator it follows that 

^(x,t) = e*^*0(x,O)e'^^* => |(/>(x);t) = e^-^V(x);0). (2.5) 

Eq. (j2.4p can then be viewed as the transition amplitude for the field to return to the same 
state after an imaginary time — when the time-development is given by the Hamiltonian 
H - inNi, 

z{T, V, fii) = [ m (0(x); t - im^y,t). (2.6) 



Dividing the time interval into infinitesimally small pieces and inserting at every point a 
complete set of position and momentum eigenstates this can be cast into a path integral 
form (for details see e.g. [34,35]) 



Z{T,V,fii) = / P^Pvrexp 



t-if3 r- 

dt' / d^x (j){yL, t')TT{x, t') - n{4>, vr) + HiMi{(f), it) 

(2.7) 

where TC and J\f are the Hamiltonian and number densities, respectively, and = dt(j)- 
When Ti. — /XjA/i is at most quadratic in canonical momenta, the momentum integration 
can be done. In gauge theory it is useful to first reintroduce the Gauss' law by treating the 
temporal gauge field component Aq as an independent variable, which, when integrated 
over, would be replaced by the stationary value satisfying Gauss' law. 

Performing the momentum integrations, we get back to the Lagrangian formulation 



Z{T,V,fii 





- i-t-il3 r 


j exp 









(2.^ 



where the integration is now over both canonical and auxiliary fields. The Lagrangian C 
usually differs from the one we started with. In particular, the momenta in Eq. ()2.7p must 
be replaced with the values solved from 

-/-(x, t) = {H[<t), vr] - /XiiVi[(/>, tt]) , (2.9) 
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so that in the end we have 



(2.10) 



Moreover, in a gauge theory one usuahy includes an additional gauge fixing term into 
the Lagrangian using Grassmannian ghost fields in order to have less constraints on the 
integration variables. 

As can be seen in Eq. (12. 6p . the partition function is computed as an integral over 
amplitudes with the same field configuration at both end points, (p{t — ip,x) = (/>(t,x). For 
fermionic variables it follows from the anticommutation properties of Grassmann variables 
that the trace has to be computed with antiperiodic condition ip{t — if3,x) = — ^(t,x) 
instead. Both boundary conditions can be verified by inspecting the two-point function, 
taking into account the correct time ordering of the fields [34] . 

When extending the time coordinate to complex values, the integration path is no longer 
unique. It can be chosen to fit the problem in question, with some minor restrictions. 
The time arguments of the operators whose thermal averages we are computing should 
obviously lie on the integration path. Also, the imaginary part of t should be nonincreasing 
in order to have a well-defined propagator. There are two conventional choices for the path, 
leading to two different ways of computing at finite temperatures. 

First, one can choose to include the whole real axis by first integrating from —to to tQ, 
then down to to — ia, with < a < /?, back to —to — icr and finally down to —to — iP, 
in the end letting to — > oo (see e.g. [36]). This approach leads to the so-called real-time 
formalism, which has the advantage that one can directly compute real-time quantities 
without having to analytically continue the final results to Minkowski space. However, 
in this formalism the number of degrees of freedom is doubled, with unphysical fields 
living on the lower horizontal part of the integration path and mixing with the physical 
ones. This in turn requires the propagators to be extended to 2 x 2 matrices, leading to 
complicated perturbation theory. We will not use the real-time formalism in this thesis. 

A simpler choice is to integrate down the vertical line t(r) = to — ir, t = 0../3, which 
leads to the so-called imaginary time formalism. The choice of to does not affect the 
results, so one can choose to = and replace the time coordinate in Eq. (j2.8p by r = it: 



The functional integral is over periodic or antiperiodic fields as described above, and the 
Euclidean Lagrangian is the same as in Eq. ()2.10p . rotated to Euclidean space with 
the replacements 



where 'E' and 'M' stand for Euclidean and Minkowski space quantities, respectively. In 
the following we will always work in Euclidean space unless otherwise mentioned, and drop 
the 'E' superscripts. In the above equation, Afj_ represents any four-vector, in particular 
the gauge field components. There is no doubling of degrees of freedom in this formal- 
ism, and for static quantities, such as the free energy or screening masses, it is usually 
simpler to compute in imaginary time. Other results have to be analytically continued to 
real time arguments, and while in principle this can be done with some mild regularity 
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assumptions [37] , in practice some additional model assumptions are required to carry out 
the continuation. 

Because the fields are required to be periodic, the imaginary time direction can be 
viewed as a closed circle with circumference /? = l/T. The momentum component in a 
compact dimension is quantized, so the fields can be decomposed in the momentum space 
as Fourier series 

, / N rr. I / \ T f "^mrT (bosons) , 

Hr, X) = T ^ ^„(x)e-- , = | ^ ^^^^ J , (2.13) 

n=— oo ^ 

where cOn are referred to as Matsubara frequencies [38]. From the gauge transformation 
rule for the gauge field components 

nA^n-^ - ^{d^,n)n-^, n{x) = exp[i5r"a"(x)] (2.14) 

it is easy to see that the gauge transformation functions a" have to be periodic as well, so 
the ghost fields will have bosonic Matsubara frequencies despite of their anticommuting 
nature. 

2.1.1 Renormalization 

The thermal environment changes the boundary conditions and the propagators from their 
zero-temperature forms. Fortunately, this does not introduce any new ultraviolet diver- 
gences, but the usual renormalization procedure remains unchanged and the counterterms 
have precisely the same values as at T = (depending on the scheme). Intuitively this is 
easy to understand, since only the excitations with wavelengths > /3 can see the periodicity 
of the time direction, while the renormalization is only concerned with divergences related 
to the short distance behavior of Green's functions. The divergence structure is then pre- 
cisely the same as in the zero-temperature theory and one can choose a T-independent 
renormalization scheme such as the MS scheme. 

To see this in some more detail, we note that the free propagator at finite tempera- 
ture can be viewed as an explicitly periodic combination of zero-temperature Euclidean 
propagators [36], 

oo 

S'ir(r,x;T)= ^ Sf{t + np,x;T = 0), 0<t</?. (2.15) 

n=— oo 

The zero-temperature ultraviolet divergences requiring renormalization arise from the 
short-distance singularities at x'^ = 0. The only term in the above sum where we can 
have = {t + n/3)^ -I- x^ = is the n = term, which does not depend on temperature. 
The divergences of the thermal propagator are therefore correctly removed by the T = 
counterterms. At higher order diagrams these divergences are multiplied by T-dependcnt 
finite parts of the diagram, so the general proof of renormalizability and T-independence 
of counterterms is somewhat more involved, but it follows from a similar decomposition 
of propagator into a singular T = part and an analytic T-dependent part [39]. 

As the parameters of the theory are renormalized, they also run with the scale according 
to the renormalization group equations. The actual equations are again the same as in 
T = theory, but the choice of renormalization point is complicated by the appearance 
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of new scales ttT and fx in addition to the external scales present in the problem, as well 
as the the scales gT, g^T generated dynamically by interactions. If these scales are very 
different, removing the large logarithms by a suitable choice of scale may prove difficult, 
and a careful analysis of the scale hierarchy is required to construct a good perturbative 
expansion. 

While the ultraviolet divergences are unaffected by the finite temperature, at the infrared 
end the situation is very different. The finite extent of the temporal direction causes the 
field components with wavelengths S> l/T to see the space effectively as three-dimensional, 
and this gives rise to many new infrared divergences. These will be treated in more detail 
in the following section. 



2.2 Dimensional reduction 

In this section we will review the rationale for dimensional reduction in the more general 
context of low-energy effective field theories. We will also discuss the finite-temperature 
infrared divergences and the resummations needed to get rid of them. 



2.2.1 Effective Lagrangians in general 

One of the fundamental properties of physics is that phenomena at some specific distance 
scale can be effectively described by a theory which does not depend on the physics at much 
shorter scales. This is fortunate, for otherwise we would not even be able to describe the 
trajectory of a thrown ball without knowledge of beyond the standard model physics. The 
same behavior, known as decoupling, is also present in quantum field theory, where it is by 
no means obvious that the heavy particles inevitably occurring as internal legs in Feynman 
diagrams can be neglected. The proof that the high-energy modes only contribute to long- 
distance phenomena by renormalization of the parameters and by corrections suppressed 
by inverse powers of the heavy masses is contained in the celebrated theorem of Appelquist 
and Carazzone [40]. Prom this point of view, every physical theory can be viewed as an 
effective theory, equivalent to the underlying more fundamental theory in some finite 
energy range. 

Formally, if the underlying theory is known, the effective theory for light modes 0; can 
be written as a path integral over the heavy modes 

giSeffM = j p<5^ expi5[0,,$;,], (2.16) 

where the effective action 5'cff[0i] is in general a non-local functional of the light fields. 
Analytically the path integral can only be computed in the Gaussian approximation around 
some given field configuration 

{^h{x) - ^h{x)) ($^2/) - ^h{y)) ■ (2.17) 



1 f^d 
+ - d'^x(ry- 



2 J S^hix)S^h{y) 



'i>h='l>h 



Choosing to be a saddle point of the action, dS[(j)i,^fi\/^^h = 0) the integration over 
$/j gives the effective action (for bosonic $/i) as 



7 



where the last term depends on cpi both directly through S and through the saddle point 
condition which makes a functional of The Gaussian approximation corresponds to 
the one-loop level in heavy-loop expansion; if we want to go beyond that the path integral 
can no longer be computed analytically, but we have to resort to perturbation theory or 
some other approximation. 

While the heavy fields can be integrated out as shown above, the resulting effective 
action is generally a complicated nonlocal functional of the light modes and cannot be 
cast in the form of an effective local Lagrangian density without some additional approx- 
imations. An often used method is the derivative expansion, where the non-local terms 
are expanded in the light field momenta p over the heavy field mass M, leading to 

5eff = / d'^xCeS + E^" (^)" ' (2.19) 

where On represent operators suppressed by powers of the heavy mass. In terms of Feyn- 
man diagrams this means that the effective action is computed with only heavy fields on 
the internal lines, since the action is made local in the light fields. The form of Eq. ()2.19p is 
precisely what should be expected based on the decoupling theorem: parameter renormal- 
izations and heavy mass suppressed operators. There is a twist, however, since the light 
particle momenta need not be small when the non-local operator is embedded in a multi- 
loop graph and interacts with heavy fields, and the derivative expansion may then fail. 
For example, in the large-mass expansion at zero temperature [39,41] it is well known that 
one needs to take into account also the diagrams with light internal lines in order to get 
the correct low-energy effective Lagrangian. This will also be the case in the dimensionally 
reduced effective theory at high temperatures, as we will show later on. 
As an illuminating example, consider a theory with two scalar fields [42], 

>c = ^a^'/'^'V - ^^202 _ + ^^,,'^>^^'<^> - ^mH^ + ^A02$2^ (2.20) 

in the limit m <^ M. This is similar to the situation at finite temperature where (p can be 
thought as the static (n = 0) Matsubara mode, while the heavy field mass is of the order 
2ttT. In this model the dependence on the heavy field is quadratic, so we can exactly 
integrate out giving 

oo , ^ 

5eff = 5[(^] + '-Tt H-d^ -M' + A02) = S[^] - ^ ^ -Tr [{d^ + M^)" V']' , (2.21) 

k=l 

where in the last step we have dropped a (/)-independent term and expanded in the small 
coupling A. The first term in the sum (k = 1) is represented by Fig. 12.1( a) and contributes 
by a local term to the mass renormalization, 

iA„ , „n , ,9s_1 ,0 iA /" ,J ,o , , /" 1 



where we have used dimensional regularization to control the ultraviolet divergence in the 
momentum integration, with the conventions 



. , , , d = 4-2e. (2.23) 

Air J J (27r)'^ ^ ^ 
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(a) ' (b) ' (c) 

Figure 2.1: Diagrams in the effective action for a tlieory witli two scalars. Solid lines represent 
the heavy field, dashed lines the light one. 



Here n is the (arbitrary) dimensional regularization scale, modified to include the constants 
typical of the MS scheme. 

The A^-term, however, already shows where the derivative expansion causes problems. 
A straightforward computation of the diagram in Fig. 12. 11 (b) gives 



X ( In — ^ — / dt m 



. (2.24) 



The first two /c-independent terms contribute to the renormalization of the 4-point vertex. 
The remaining logarithm is a non-local operator connecting two 0^ products at different 
points. For small k'^ the integrand can be expanded in k'^/M'^, leading to a series of local 
four-point derivative couplings of the form (/)(9^/M^)"(?i>. However, when this operator is 
part of a larger diagram there is no guarantee that k^ is small. 

For example, the diagram in Fig. 12.1( c) with one light and two heavy internal lines 
is not produced by the effective theory expanded this way. All loop momenta can be 
large, and therefore the expansion in A;^/M^ is not reliable. Computing this diagram is 
rather nontrivial [43] , but one can show that if the ultraviolet divergences are removed in 
the MS scheme, the diagram does not vanish in the limit M — > oo. To have an explicit 
decoupling where all graphs containing heavy internal lines are suppressed one should use 
a renormalization scheme where the counterterm is the negative of the graph expanded in 
the light masses and momenta [39]. At finite temperatures this may be difficult because 
of the additional infrared divergences. Moreover, we would prefer to use the MS scheme 
where the counterterms are already known to high order and have a simple structure. 

Because of the difficulties in integrating out the heavy fields as described above, at higher 
orders it is usually safer to construct the effective Lagrangian explicitly by matching the 
Green's functions. The decoupling theorem states that in a renormalizable theory the 
parameters in the effective theory can be chosen in such way that the Green's functions 
of light fields differ from those computed in the full theory by terms suppressed by powers 
of the heavy mass, 

Gn{pi,-- ■ ,PN]9,G,m,M,fi) = {0\T(p{pi) . . . (I){pn)\0) lull 
= z-^'/'GUpi, ...,PN;9*,m*,fi)[l + (l/M^)] 
= z-^/2(O|T0*(pi) . . . </.*(p7v)|0)efr [l + O (l/iVr )] , (2.25) 

where M and G are the masses and couplings in terms involving heavy fields, while 
those for terms with only light fields are labeled m, g. The corresponding effective theory 
parameters are m*,g* and 4>* = z^/'^cj). We can use this information directly and write 
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down the most general light mode Lagrangian which respects the symmetries of the original 
theory, and then compute a number of A^-point functions (usually = 2, 3, 4 is enough) 
in both theories at some conveniently chosen external momenta to fix the parameters. We 
will see more detailed examples of this procedure in the following section. 

2.2.2 Three-dimensional effective theory at high temperature 

Field theories at finite temperature contain many new mass scales in addition to those 
given by the parameters of the zero-temperature Lagrangian. Besides the temperature 
itself there are dynamically generated scales related to collective modes and screening 
phenomena, and the particle masses are modified by thermal effects as well. Renormal- 
izing the theory in the minimal subtraction scheme gives rise to logarithms of the type 
ln(m^/^^), where m can be any of the different scales in the theory. In particular, large 
scales do not decouple but instead give contributions that grow logarithmically with the 
scale. This seems to make perturbation theory useless in theories with vastly different 
mass scales, since we cannot choose a renormalization scale that simultaneously makes all 
the logarithms small. As a result, terms in the perturbative expansion contain powers of 
large logarithms in addition to small coupling and need not decrease at higher orders. 

To be more specific, in gauge theory the electric and magnetic screening scales are of 
order gT and g^T, respectively, and thus there is a clear hierarchy of scales in the small 
coupling region where we would want to use perturbation theory. The solution is, as 
discussed above, either to use a more complicated renormalization scheme or to formulate 
an effective theory and continue using the MS scheme [44]. As it turns out, it is simpler 
to carry out the computations using the effective theory. We will mostly concentrate on 
gauge theories in what follows, in particular on QCD and electroweak theory. 

In the imaginary time formalism we can write the four-dimensional theory in terms of 
the Matsubara modes of Eq. ()2.13p . For generic bosonic and fermionic fields the free part 
of the action (without any chemical potentials, although they could easily be included) is 

oo „ 

= r ^ d^x4>i {-df + [(27rnr)2 + ml]) 0„ + i^n[i{2n + l)7rr7o + ^^^^ + m/lV'n , 
n=-oo (2.26) 

which can be viewed as a three-dimensional Euclidean theory of an infinite set of fields 
with masses = + w?. If the temperature is much higher than the particle masses, 
we can use the arguments of the previous section and try to formulate an effective theory 
for the light modes with M„ ^ T, or the bosonic zero-modes since they are the only 
modes with = 0. This theory loses all dependence on the (imaginary) time coordinate, 
so we have effectively reduced the number of dimensions to three. From the point of view 
of modes with wavelengths much larger than 1/T the finite temporal direction of length 
(5 has shrunk to a point. 

While the dimensionally reduced theory cannot give any information about the time 
dependence of the theory, for static Green's functions the effective theory gives correct 
results up to corrections of order m^/(7rT)^, where m is any of the light masses. Note 
in particular that at high enough temperatures the highest unintegrated mass scales are 
the dynamically generated scales ~ gT, so the corrections to the effective theory are 
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comparable with the higher orders of perturbation theory and both have to be taken into 
account to get a consistent perturbative expansion. To gain control over which operators 
to include, power counting rules have to be established for given momentum region. At 
higher orders it will be necessary to include nonrenormalizable operators into the effective 
theory, especially if one wishes to have a theory that produces all static Green's functions 
to given order. In many cases, like when computing the free energy, it is sufficient to use 
only the couplings present already in the original theory, in which case the effective theory 
is super-renormalizable because of the lower dimensionality. 

The main advantage in using an effective theory at high temperatures is in the infrared 
physics. In general, if the theory contains masslcss bosonic fields one expects more severe 
infrared singularities when going to finite temperature, since the Bose-Einstein factor in 
real-time propagators behaves as 

This can be also understood in the imaginary time formalism, where the zero MatsTib- 
ara mode behaves like a massless particle in three dimensions, and lower dimensionality 
generally makes the infrared behavior worse. It is well known that in Yang-Mills theories 
perturbation theory at finite temperatures suffers from many infrared problems, becoming 
finally completely non-perturbative at 0{g^) [45,46]. These problems are related to mass- 
less particles, in particular to the gauge fields, whose screening by medium effects is not 
correctly reproduced by the naive perturbation theory. By definition, the dimensionally 
reduced theory has the same infrared limit as the original theory, while being computa- 
tionally simpler. The leading order contribution coming from scales of order T can be 
included in the parameters of the effective theory via the matching procedure, which is 
infrared safe, and the infrared peculiarities can then be studied in a simpler setting. In 
particular, the dimensionally reduced effective theory does not contain any fermionic fields, 
which makes it easier to study non-perturbatively using lattice simulations. 

The electric screening effects can be included by reorganizing the perturbative expan- 
sion. Computing the one- loop self-energy of a zero-mode gauge field component A^, we 
find that in the limit of vanishing momentum it behaves as 

^nuioon = 0, k ^ 0) a g'^T^S^oSuO ■ (2.28) 

The temporal component develops a thermal mass of order gT, while the other components 
remain massless. In the soft limit where k < gT it is not consistent to treat this self-energy 
as perturbation, but it should be included in the propagator instead. This means that we 
should sum all diagrams with an arbitrary number of self-energy insertions on the temporal 
gluon line to get consistent 0{g^) results, which is often referred to as resummation. In four 
dimensions one has to be careful not count any diagram twice because of this summation; 
usually this is done by adding and subtracting a term containing the self-energy in the 
Lagrangian, 

C = Cq + Ci = {Cq + 5C) + (£/ - 6C) (2.29) 

and treating the subtracted term as an interaction. In the dimensionally reduced theory 
the resummation is simpler, since the thermal mass for Aq comes out naturally from the 
matching procedure. Moreover, there is no risk of double counting diagrams, since the 
thermal mass is only created by n 7^ and fermionic modes (the mass can be computed 
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in the k = hmit, and the dimensionless graphs vanish in dimensional regularization) , 
which are not present in the effective theory. Note that the electric mass does not break 
the remaining gauge invariance, since when restricting to bosonic zero modes only we are 
also forced to only consider r-independent gauge transformations. The transformation 
rule in Eq. (j2.14p then boils down to 

= 2 TV AqT"" ^ 2 Tr nAQn^^T"" = exp[ic/a'=<J , (2.30) 

so in the three-dimensional theory becomes a massive scalar transforming in the adjoint 
representation of the gauge group. The remaining gauge invariance in three dimensions 
prevents the spatial gauge field components from developing a mass term. 

In the magnetic sector there are infinitely many diagrams that all contribute at order 
g^, and, unlike for the electric mass, they appear with so different and complex topologies 
that they cannot be resummed in a simple way to tame the infrared singularities. In fact, 
there is no gauge-invariant magnetic mass term that could be included in the Lagrangian 
for perturbatively computing beyond 0{g^), but instead the magnetic screening has to 
be treated non-perturbatively. In the very low momentum region the fields with thermal 
masses ~ gT can be integrated out as well, leaving a three-dimensional pure gauge theory 
with coupling g^ = g^T, which is the only dimensionful parameter in the Lagrangian. In 
this theory there is no small dimensionless parameter to do perturbation theory with, but 
the infrared dynamics of nonabelian gauge theory is inherently nonperturbative. 

To see how the matching of parameters in the dimensionally reduced theory goes in 
practice, we will take a closer look at the mass parameters, following to some extent [11,23]. 
The masses can be found by comparing the static two-point functions computed in both 
theories. For simplicity, we will use a scalar particle with a small zero-temperature mass 
m < gT as an example and work to order g'^, which is sufficient for many computations, 
in particular for determining the free energy to order g^ as in [2,3]. 

In the full theory the inverse propagator can be written as 

k'^ + m^ + n(A;2) = k^ +m^ + U{k^) + Uo{k^), (2.31) 

where n(/c^) includes the diagrams with at least one heavy internal line, while no(A;^) is 
the contribution of n = modes only. In the effective theory the same function reads 

k'^ + ml + U3{k^). (2.32) 

The contribution coming from the non-static modes, n(A;^), is of order g'^T'^, and the 
matching has to carried out in the region where the effective theory is valid, k < gT. Since 
integration over massive modes is infrared safe, the renormalized self-energy n(/c^) has no 
infrared divergence and can be expanded in k'^ /T'^, 

n(fc2) = n(0) + k'^m + O (^g'^^ , (2.33) 

where the terms left out are of order g^T^. Further expanding each term in loop expansion 
with coupling g, 

oo 

TT(A:) = ^tf"^(A;), where TT^"^(A;) ~ 0(^2"), (2.34) 

n=l 
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the inverse propagator in Eq. (|2.3ip reads, including terms up to 0{g'^), 



^'^^^ dF"^'^^°^ ^ + + n(^>) +rf')(0) + Hoik'). 



(2.35) 



The massive modes correspond to poles in the propagator, or the zeros of the inverse 
propagator, so we set the expressions in Eqs. (I2.32|) . ()2.35p equal to zero and solve for k'^. 
Equating the pole locations in both theories, we find the matching condition 



ml + Us{k^ 



1 



dk^ 



m^ + u^^\o) + n^^\o) + Uoik^) 



(2.36) 



By construction, the infrared behavior contained in the soft self-energies Ho and is the 
same in both theories, so this relation is infrared safe. The difference is of order g^, 



n3(A;2) = no(fc2) [1 + Oik^/T^ 



Us{k 



2\ 2 



(2.37) 



so, working at order g^, we can drop all terms containing 110,113 from the matching 
condition. We are then left with an equation for the three-dimensional mass parameter 



m 



2 - m2 + n^^) (0) + (0) -(m^ + n^') (0) ) -^n^^^ (o) . 



dA:2 



(2.38) 



As a by-product we also found the field normalization factor to order , since from looking 
at the coefficients of fc^ in both propagators we can write 



1 
T 



(2.39) 



The factor 1/T here stems from the overall factor T in Eq. (|2.26p . which is conventionally 
absorbed into the fields and couplings of the 3d theory. 

It should be noted that Eq. ()2.38p only contains contributions from the heavy scale 
r, whereas the infrared sensitive parts Ho and Us drop out. The mass parameter ms 
regulates the infrared behavior of the dimensionally reduced theory, but it is a completely 
perturbative quantity and should not be confused with the actual screening lengths that 
are sensitive to infrared physics. In particular, the thermal mass of the adjoint scalar 
in the dimensionally reduced theory agrees with the electric screening mass m^i only at 
order g'^, beyond which mei becomes sensitive to the magnetic screening [47], while on 
its part is given to 0{g'^) by the completely perturbative expression in Eq. (|2.38p . 

Apart from the gauge fields, the only other elementary boson in the standard model 
is the Higgs field, which has a negative mass parameter — z/^ in the phase of unbroken 
SU{2) X U{1) symmetry. Near the electroweak phase transition the Higgs field mass is 
a special case in the power counting, since the T = mass parameter and the thermal 
corrections almost cancel each other, giving 



-1/2 + g'^T'^ 



3rp2 



g'T 



or smaller, depending on how close to the phase transition we choose to work. To have 
a better separation of scales, it is necessary to integrate out the fields with masses ~ gT 
when computing close to the electroweak phase transition, as we did in [3]. This leads to 
a theory containing only the Higgs field and spatial gluons. The thermal mass m|(T) is 
the leading term in the Higgs field effective potential, which drives the phase transition. 
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The above matching computation gives another example of how the expansion in loops 
and momenta can be identified when the the correct momentum region is known. At high 
temperatures, the mass parameters can be estimated as gT and the momenta at most of 
the same magnitude, in the region where dimensional reduction is valid. The required 
level of matching is determined by the problem in question and the accuracy goal one 
wants to reach. For example, for computing the free energy to order we needed the 
couplings only at tree-level, but the mass parameters to two-loop {g'^) order. 

A more general analysis given in [23] states that in order to have a theory which gives 
the same light field Green's functions as the full theory up to corrections of order 0{g'^), 
we need to match the parameters at least to this order. To be more precise, the coupling 
constants are required to one-loop level g^ = T[g^ + g'^) and adjoint scalar (temporal 
gauge field component) masses to two-loop accuracy m| = T^{g^ + g*). If the theory 
contains a light scalar field such as the Higgs field, its thermal mass should be computed 
to three-loop level rrig = —v'^ + T^{g^ + 5^ + 9^), since the first terms cancel each other, 
and the mass is of order g'^T^ close to the phase transition. The same analysis shows that 
beyond 0(5^) it is necessary to include non-renormalizable 6-dimensional operators into 
the effective theory. 

Apart from the simple power counting, the importance of the higher order operators 
inevitably resulting from the reduction step is difficult to estimate. In [23, 48] it is noted 
that in both abelian and SU(2) Higgs models these operators are further suppressed by 
small numerical coefficients in addition to powers of the coupling constant, and thus give 
only very small contributions. The operators following from the second reduction step, 
where the scales ~ gT are integrated out to give a pure gauge theory, can be consistently 
treated as perturbations with respect to the tree-level Lagrangian, as discussed in [49]. 

For matching purposes we still need to compute some Green's functions in the full theory, 
but using the effective theory this only has to be done once, after which the computations 
can be carried out in the simpler effective theory. For both QCD [9, 11] and electroweak 
theory [23] the matching has been carried out explicitly to order g^, and for a generic 
theory containing scalars, fermions and gauge fields the rules given in [23] can be used to 
find the parameters of the effective theory. The QCD coupling has even been matched to 
two-loop [gl = Tig"^ + g"^ + 5^)] level in [50]. 

While the effective theory approach saves us from computing multiple complicated sum- 
integrals, at finite temperatures the main advantages of dimensional reduction lie in the 
easy way to organize the resummations and separating the contributions of different scales. 
Eventually non-perturbativc methods such as lattice simulations are needed to handle the 
infrared limit, but the dimensional reduction methods allow us to work out the parame- 
ters with completely perturbative methods, and then apply the computationally intensive 
methods to the simpler three-dimensional theory. Lattice simulations in the dimensionally 
reduced theory are easier because there is one spatial dimension less, no fermions and the 
shortest scales < l/T have been integrated out. 
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Chapter 3 

Pressure of the standard model 



At high temperatures the local SU{2)l x U{1)y gauge symmetry of electroweak theory 
is restored. The phase transition is driven by the Higgs field, whose effective potential 
is modified by thermal corrections in such way that the vacuum expectation value of the 
field vanishes when the temperature is raised. Because of the possibility of the phase tran- 
sition being strongly first order and contributing to the baryon number asymmetry, the 
effective potential has been extensively studied both by 1-loop [51-53] and 2-loop [54-56] 
perturbative calculations and by dimensional reduction [23] combined with lattice simu- 
lations [24-28]. In those works it was shown that in the standard model the electroweak 
phase transition is a crossover for realistic Higgs masses. 

Apart from the effective potential computations, the thermodynamics of electroweak 
theory has not been studied in detail. In [2, 3] we computed the most fundamental ther- 
modynamic quantity, the free energy, for electroweak theory at high temperatures. This 
computation is very similar to the evaluation of the free energy in QCD, with the main 
differences coming from the presence of a light scalar field driving the phase transition 
and the multitude of scales and couplings leading to a very complicated general struc- 
ture. Together with the QCD result and the few terms mixing the strong and electroweak 
couplings, this computation gives us the free energy of the full standard model. 

Partial derivatives of the free energy give the basic thermodynamical quantities as in 
Eq. ()2.3p . It should be noted here that we are computing in the grand canonical ensemble, 
whose partition function gives the grand potential 0, = — TlnZ, but at zero chemical 
potentials this can be identified with the free energy F = Q+^iNi. In the thermodynamical 
limit y — > oo the free energy density equals the pressure, F = —pV, so for simplicity we 
will we talking about pressure from now on. 

The energy density and pressure are particularly interesting, since they control the ex- 
pansion of the universe at its very early stages. Temperatures higher than the electroweak 
crossover cannot be reached experimentally, but they were present in the early universe. 
The relic densities of particles decoupling from the ordinary matter are sensitive to the 
evolution of the universe, which in turn is governed by the equation of state. Recent 
measurements of the cosmic microwave background suggest a sizeable amount of cold 
dark matter, which could be explained by weakly interacting massive particles (WIMPs) 
(see [57] for a review). Given a theory describing WIMPs, we need to know the evolution 
of the universe at the time of their decoupling as well as at later times to make predictions 
of the present situation. In [58] it is estimated that a 10% change in the equation of state 
leads to 1% difference in relic densities, which is visible in future microwave observations. 



15 



3.1 Perturbative evaluation of the pressure 



The electroweak sector of the standard model is given by the Euchdean Lagrangian 



+ QLlPqi + URipUR + dRipdR + igv qLT^^*tR - tR{^*yT^qL 



(3.1) 



where GJJ^ 



duAl + ge 



abc Ab 



A^Al and F^^ 



duBy — dpBn are the field strengths of 



the weak and hypercharge interactions, $ is the Higgs field and the covariant derivatives 
act on the chiral fermion fields and the Higgs field as usual (for details, see Eq. (2.3) 
of [2]). We only include the Yukawa coupling for the top quark, since for other particles 
the Yukawa couplings (which are proportional to particle masses in the broken symmetry 
phase) are orders of magnitude smaller. 

When the Euclidean action is given, the pressure can be computed as the logarithm of 
the partition function. 



piT) 



T 

lim — In 



VAVi/jVipV^ exp 



/3 



dr d^xC{A,^l;,^,<^) 



(3.2) 



where the path integral is over all fields in the Lagrangian. As described in the previous 
chapter, a straightforward perturbative evaluation of the path integral Eq. (j3.2p fails 
because of infrared divergences. The solution is to resum a class of diagrams by means of 
an effective theory, using dimensional reduction. 

In the first level of dimensional reduction all non-static modes, in particular all fermions, 
are integrated out. This leads to an effective theory S'e, whose parameters are matched 
by perturbative computations in the full theory with no resummations. 



PiT) = pe{T) + lim ^ In / VAkVAoV<l> exp (Se) ■ 

V^oo V J 



(3.3) 



Note in particular the appearance of parameter pe{T), which is the contribution of the 
non-static modes, or scales ~ vrT, to the pressure. This parameter can be also viewed as 
the matching coefficient of the unit operator by looking at the (unnormalized) expectation 
value of unit operator in both theories, 



(1) 



full 



Tr 1 • pf, 



ull 



full 



-F/T 



-Fe/T+IuZe 



(3.4) 



where F = —pV. Since the matching is infrared safe, all parameters of S'e and also 
Pe are series in g^. In addition to curing some of the infrared problems, this approach 
makes full use of the scale hierarchy T ^ gT ^ g^T by separating the contribution from 
each scale into successive effective theories, whose contributions enter at different levels of 
perturbation theory. For example, it is easy to see that the dimensionally reduced theory 
S'e in Eq. (13. 3p starts to contribute at level Tm| ~ g^T"^. 

The effective theory Se still contains two different scales gT and g'^T, the latter of which 
is related to non-perturbative magnetic screening effects. If one wishes to go further using 
perturbation theory, it is useful to integrate out the electric scales gT as well, giving 



p{T) = pe{T) + pm{T) + f VAk exp (-Sm) 



(3.5) 
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where Sm only contains the spatial gauge fields. Close to the phase transition this step is 
more complicated because the scalar mass is very light, and deserves a separate discussion 
in section 13. 2[ The only dimensional parameter in 5m is the gauge coupling « g^T, 
so this theory begins to contribute at order Tg^ ~ g^T^^ this term being completely non- 
perturbative. We have only kept terms of order cf' in our calculation of the pressure, so 
this non-perturbative contribution can be dropped. The purpose of this second reduction 
step is that the computation of pM in the first effective theory Se can be considered as 
a matching computation, without having to worry about the resummations needed for 
spatial gauge fields. 

The electroweak theory contains many dimensionless coupling constants, and we need to 
establish a power counting between them in order to determine which terms to include in 
the perturbative expansion. We have decided to use the weak gauge coupling as reference, 
and, denoting the strong and hypercharge couplings by g^ and 5', respectively, make the 
simple choice 

\ /2 2 2 2 2 ^ 2rn2 In a\ 

A'-S' ^ Qs ^ 9y ^ 9 , ^ ^9 T , (3.6) 

which corresponds to three-loop expansion in all couplings. Numerically this is not the best 
choice, since the strong and Yukawa couplings are large compared to electroweak couplings, 
and we underestimate their importance. However, trying to incorporate higher orders of 
gY would require four-loop sum- integrals, which we do not know how to perform. The 
strong coupling is even harder, since the g^ order suffers from the same infrared problems 
as any nonabelian gauge theory. In [23] the rule g''^ ~ 5^ is used, but there is no danger 
and practically no extra work in keeping terms of order g'^ as well. It should be also kept 
in mind that the one-loop renormalization group running of the couplings is such that 
g'(T) grows with temperature, whereas g{T) decreases. 

After all the preparations are done, it remains to actually compute the pressure. We 
start by evaluating p-^ to three-loop order. At the 4d full theory level the Higgs mass 
parameter i/^ is treated as a perturbation, so we expand the propagators in z/^. This is 
possible since the matching procedure is infrared safe. The resulting massless sum-integrals 
can be evaluated using the methods developed in [9] and can be conveniently read from 
the Appendix A of [11]. The largest work lies in writing down all the required diagrams 
with correct symmetry and group theory factors, reducing them to integrals given in [11] 
and summing everything together. Note that the diagrams with only static modes do not 
have to explicitly subtracted, since they vanish in the dimensional regularization due to 
lack of dimensionful parameter. 



Schematically, the generic form of pE is 
Pe{T) = T 



V 



4 



+ ^«i?.. + r^-0(5'), (3.7) 



where the summation is over gf = g^ ^ g'^ ^gl^gy,^ and the values of all nonzero coefficients 
(all combinations except aEXsidEsu, and aEsiOiEss since we exclude the pure QCD terms) 
can be found in the Appendix A of [2]. To keep track of different contributions they are 
given in terms of the group theory constants, which for SU(2) read Tp = 1/2, C-p = 3/4, 
d-p = 2, Ca = 2 and = 3. 
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We have normalized the pressure so that (the real part of) the pressure at the symmetric 
phase vanishes at zero temperature, p{T = 0) = 0, in order to exclude the large vacuum 
energy contribution and the related divergences. This normalization is already taken into 
account in Eq. (13. 7p . where we have subtracted a term proportional to i/^ computed at 
zero temperature. The T = computation differs from the high-temperature expansion 
in Eq. ()3.7p . and the difference is contained in the remaining coefficient aEuu- 

The renormalization of the 4d parameters does not remove all divergences, but 1/e 
terms can be found in most of the coefficients of 0{g^) terms, corresponding to infrared 
divergences that cancel against similar terms in pM- Having stated above that the matching 
computation is infrared safe, we should elaborate on the nature of these divergences and 
their cancellation a bit more. Dimensional regularization simultaneously handles both the 
infrared and ultraviolet limit, and it is not easy to tell the divergences apart. 

The electroweak theory is known to be renormalizable. The computation of pE cannot 
thus contain any ultraviolet divergences, since they are removed by the counterterms. 
However, there are diagrams that are both ultraviolet and infrared divergent and vanish 
in dimensional regularization. We can see how they behave through the following simple 
example. Consider the logarithmically divergent integral 

k^~~\ Att J J (27r)4-2^A;4' 

which vanishes in dimensional regularization. This can be written as a sum of two integrals, 
one divergent at the ultraviolet and the other at the infrared momenta, 

1 m2 f , 1 2m2(l -x) 



fc2(/c2 + m^) k^(k^+m'^) JkJo [k"^ + xm?]^ [k"^ + xm? 
1 /eV 



I f I 1 



^ dxT{e)x-^ + T{l + e){l- x)x- 



167r2 



(3.9) 



The renormalization counterterms remove the ultraviolet divergence 1/euv here, leaving 
the infrared divergent part — l/eiR. The vanishing diagram therefore contributes with an 
infrared divergence when renormalized. 

The scaleless diagrams at finite temperature are precisely those with only static modes, 
and no summations over the Matsubara frequencies. The dimensionally reduced effective 
theory contains the same diagrams, but with self-energy corrections resummed to give 
masses on some propagators. These masses do not change the ultraviolet behavior of the 
diagram, but regularize the infrared limit, so the divergence structure is just l/euvi which 
precisely cancels against — 1/eiR, from the full theory computation. 

The contribution of the scales gT can be calculated from the path integral in Eq. ()3.3p 
once the dimensionally reduced theory is known. Before matching, we need to consider 
the most general renormalizable (in 3d) Lagrangian respecting the symmetries of the full 
theory, 

Se = j d?x \g%G% + -^FijFi, + (A1>)^(A^) + mi<Dtci> + X^{<^^^f 

+ \{D,Alf + \mlAlAl + \\A{AlAlf + \{d,B,f + \m%B,B, + \\bBI 
+ h^^^^AlAl + h'^^^^B^B^ - ^5353^o^^^Sr"$ , (3.10) 
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where we have included masses and quartic self-interactions for the scalar fields Aq,Bo, 
the former temporal components of the gauge field. The field Aq transforms in the adjoint 
representation of SU(2), wheres Bq does not interact with the gauge fields due to the 
abelian nature of U(l). 

In order to relate the pressure pM to the full theory we need to know the parameters of 
the effective theory in terms of the full theory couplings. For the couplings the leading 
order results are sufficient, since the two- loop diagrams in the effective theory are already 
of order (7^. The matching then boils down to absorbing the factor of T^/^ to couplings to 
give them the correct dimensions, 



The quartic couplings for the adjoint scalars are of higher order than we need in our 
computation. 

The matching of the mass parameters is more complicated. Since the leading order 
(one-loop) diagrams are of the order Tm^ ~ g^T^, we need order g^ terms in the expres- 
sions for the masses to get the pressure up to 0{g^). Moreover, the two-loop diagrams 
contain ultraviolet divergences, so we need also the 0{e) terms for the masses when using 
dimensional regularization. The mass parameters in electroweak theory have already been 
computed in [23] apart from the g'^e terms, which we have evaluated in [2]. The masses 
are found by matching the two-point functions at vanishing external momentum as in 



and similarly for m^. The coefficients f3Ex can be found in the Appendix B.l of [2]. It 
should be noted that there are no divergences in these coefficients, but the renormalization 
of the 4d theory is enough to make the adjoint scalar masses finite. In the 3d theory these 
parameters are renormalization group invariant to this order, and only start running at 
order g^ with terms proportional to X\ and g^XA [55]. 

The fundamental scalar mass has 1/e divergences that are not removed by the renor- 
malization of the full theory. These are again related to the infrared limit of the static 
modes and are removed by the counterterms in the effective theory. Looking from the 
dimensionally reduced theory, the matching procedure produces the bare mass which we 
can either split into the renormalized mass and counterterms, or continue using the mass 
parameter with 1/e terms included as we did in [2] . Either way, the divergences will cancel 
in the final result for the pressure, and itself is not a physical parameter we would be 
interested to study in detail. 

Using the MS scheme in the effective theory to renormalize the Higgs mass, we get the 
finite result 





Eq. (lOHD . 

The general form of the adjoint scalar masses is 




(/3ij3 + 0(e)) + 0(/) 




mi(A) 




) 
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PAA + .r, PBB + -n~T9PAB + PAX + T1~^PB\ + no PAA 



(47r 



(4vr) 



(4vr) 



(4vr) 



(4vr) 



(47r)2 



/2 2 2 2 

9 gy o , 9s9y a , 



Ag^ 

(47r)2^^^ ' (47r)2'^'''' ' (4^)2^"' ' (47r)2'"^^ ' (47r^2 



/?Ay + 



5^ 



(3.13) 



which depends on the MS renormahzation scale A replacing /u in Eq. p. 230 . The param- 
eters are linear combinations of C,'{—1), je and In(A/47rr), and they are given explicitly 
in the Appendix B.2 of [2]. The mass counterterm can be read from the matching com- 
putations, 



5mo 



rp2 



A , /4 , 1^ 2 /2 

+64^ +32^^ 



Ixg' - Ixg" + 3A2 



(47r)2e 

This coincides with the counterterm computed directly from the effective theory 5e 



(3.14) 



1 



(4^)2e 



39 ^ ^ /2 ^^2/2 

64^'^ ~'~ 64^'^ 32^^^^ 



^Asgi - ^Xsg'i + 3Ai + hi - Shsgl + 2h'i 



(3.15) 

when the relations ()3.1ip between couplings are taken into account. Since the 3d theory 
is super-renormalizable and has only a finite number of divergent graphs, the counterterm 
in Eq. p. 150 is actually an exact result, without any higher order corrections [55]. 

Since we are computing only vacuum diagrams in the effective theory, matching the 
fields is not required, apart from what was included in the mass parameter computations. 
All the required parameters are then known, and the pressure pM can be computed by 
evaluating all one-particle irreducible vacuum diagrams up to three-loop level in terms of 
these parameters. Apart from a gauge boson loop with both Higgs and Aq self-energy 
corrections, all the required 3d integrals are computed in [11]. Because of the massive 
propagators the general structure of the result is much more complicated than for p-^: 



PujT) ^ _1_ 
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2n3/2 
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in 
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(3.16) 
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The coefficients B^yz and the coefficient function b{x) are hnear combinations of 1/e, In 2, 
TT^ and ln(//3/M), where /^a is the 3d dimensional regularization scale and M can be any 
combination of the different mass parameters mo , mj-, , . The detailed expressions can 
be found in the Appendix C of [2]. 

The part of the pressure coming from the electric scales in Eq. (j3.16p has many new 
features that are not present in the corresponding computation for QCD. In particular, 
the only dimensional parameters in dimensionally reduced QCD (known as EQCD) are §3 
and rriD, so the possible terms are, for dimensional reasons, rrijy, g^mi) and g^mY), with 
divergent coefficients containing ln(/i3/mD)- This is in sharp contrast with the abundance 
of different terms in Eq. (|3.16|) : not only are there many combinations of couplings and 
masses, but also completely new kinds of expressions like mf/rrij and ln{mi/mj). 

All the coefficients Bxyz and b{x) of the 0{g^) terms have ultraviolet divergences, but 
they cancel against the 1/e terms in the mass counterterm (Jmg . The renormalized mass in 
Eq. (|3.13p depends on the renormalization scale A through logarithms ln(A/47rT), which 
come with the divergences as usual. They cancel against the corresponding logarithms 
ln(//3/M) in Eq. (j3.16p . leaving terms like g^m-£i\n.{m-£i /T) ~ 5^ In g. This kind of terms 
are not present in EQCD, where all the mass parameters are finite at 0{g^). If we choose 
= A, the scale dependence in pM vanishes completely at 0{g^) and 0{g^) when the 
running of the 4d couplings is taken into account. 

The remaining 1/e terms shown explicitly at two- loop level in Eq. (j3.16p cancel against 
the infrared divergences in p-^, Eq. (j3.7p . Also there the cancellation between terms coming 
from scales vrT and gT results in large logarithms of order 5^ In g. These terms are also 
present in the QCD pressure, where they were originally derived from the requirement 
that the pressure should not depend on the scale at 0{g'^) [8]. 

The presence of terms proportional to ln{g) shows that we cannot choose the scale 
in such way that the large logarithms would completely vanish. The use of an effective 
theory to separate the contributions from different scales is often advocated by the absence 
of large logarithms, but as we see, the infrared divergences mix the different scales, and 
logarithms of ln{gT/T) are left in the final result. A stronger argument for formulating the 
problem in terms of effective theories is the proper handling of resummations and isolating 
the non-perturbative infrared behavior into a simpler theory, as discussed in section [2.2.21 

3.2 Pressure near the phase transition 

One of the most interesting properties of the electroweak theory is the crossover phase 
transition. When the temperature is lowered, the effective potential of the Higgs field 
develops a new minimum at some finite value (<I>) 7^ 0. This phase transition gives masses 
to all quarks and leptons (except neutrinos) as well as to those gauge bosons mediating the 
weak interactions that correspond to the broken part of the eletroweak SU(2)i x U(l)y 
symmetry. 

In the deconfinement phase transition of QCD the strong coupling grows very large 
and perturbation theory cannot be used to study the phase transition. This is not the 
case in the electroweak symmetry breaking, since the confinement radius of weak interac- 
tions is tremendously large {g has a Landau pole at 1/Aew ~ 10^ m) compared to the 
relevant distance scales at the electroweak transition temperatures (l/Tc ~ lO"^*^ m for 
Tc ~ 200 GeV). We can then try to use perturbation theory in studying the electroweak 
phenomena close to the phase transition, in particular to extend our previous computation 
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of the pressure down to transition temperatures. Besides the pressure playing a central 
role in describing thermodynamics near the transition, this also allows us to test the va- 
lidity of dimensional reduction at the phase transition, something that cannot be done in 
QCD. 

Note that perturbation theory is not able to describe the transition itself correctly with 
physical Higgs masses. Both two-loop effective potential calculations [55], although limited 
to ruH ^ iTT-w, and the e-expansion analysis extrapolated to e = 1 [59] suggest the presence 
of a first order phase transition for large Higgs masses, while in reality the transition is 
of crossover type for mn ^ 72 GeV [25,27,28]. In the unbroken phase the Higgs field 
has a non-vanishing vacuum expectation value, which we would have to incorporate in 
our computations to build the perturbative expansion around the true, physical vacuum. 
Neverthless, because the coupling constant stays small in the transition we are able to 
perturbatively compute physical quantities while approching the transition from above, 
as long as we stay in the symmetric vacuum. 

The pressure computed in Eqs. ()3.7p . (l3.16p cannot be directly continued down to the 
phase transition, because it becomes singular as we approach the critical temperature 
(which exists in perturbation theory). The reason for this is that the thermally corrected 
Higgs mass in Eq. (|3.13p becomes very small at the phase transition, and finally turns 
negative a little below the transition. At leading order it is this negative mass parameter 
that makes the symmetric phase unstable and drives the phase transition. The singular 
effects of small thermal Higgs mass can be seen in pM, Eq. (|3.16p . which contains terms 
like m^/ms and mo ln(/Li3/m3). 

The origin of this problem is in our power counting, which assumed ~ mj) ~ gT. 
Close to Tc this assumption fails and we run into the same infrared problems as in the 
original theory, since the resummation of non-static modes no longer gives a finite mass to 
the Higgs field. In Fig. 13. II we have plotted the ratio of the renormalized Higgs mass to the 
SU(2) adjoint scalar mass in both the full standard model (including strong interactions) 
and the SU(2) + Higgs theory which we have used as a weakly coupled toy model. The 
numerical values of the parameters are given in section 13.31 and the renormalization scale 
for 7713 is chosen as A = 2ttT. As the figure shows, the mass ratio drops steeply as the 
temperature approaches the phase transition, and at T — Tc < 50 GeV the assertion 
m3 ~ mD clearly fails. 

The solution is to resum yet another class of diagrams, the adjoint scalar self-energy 
corrections on the fundamental scalar line. We then have one more level of dimensional 
reduction, and instead of Eq. ()3.5p the pressure is given by 

p(T) = pe(T) + PMi(T) + pu2{T) + f In y VAkV<^ exp {-Sm) , (3.17) 

where pE is the same as before, Eq. (j3.7p . and pMi is computed from the effective theory 
5e in Eq. (I3.10p . treating now the light scalar mass as a perturbation and expanding 
the integrals in mg/m^. The contribution from scales rn^ < g'^T^ is contained in the con- 
tribution pm2 ) which is computed from an effective theory containing only the fundamental 
scalar field and the spatial components of the gauge bosons, 

Se2 = j d'^x \g%G% + -^FijFi, + (A^)^(A^) + mi$^$ + M^^^? , (3.18) 

where the gauge couplings and the scalar self-coupling do not get any matching corrections 
at this level, (^1, g'^, A3) = (^f^, g'^, A3). The mass parameter, however, now also resums 
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Figure 3.1: The ratio of the fundamental scalar thermal mass to the SU(2) adjoint scalar mass. 
The regularization scale is chosen as A = 27rT. 



the adjoint scalar loops, which we have to compute up to one-loop level {g^ in our power 
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Apart from the 0{e) terms, this expression has been previously computed in [23]. The 
correction is of the same order as the leading term m^, showing that resummation is 
necessary to get consistent results. 

For simplicity, we have neglected the order g'^ corrections to the scalar mass. This can be 
justified by power counting arguments, since for m| ~ g^T^ the two-loop corrections would 
contribute parametrically at order m^gl ~ g^^/'^T^, which is strictly speaking higher than 
0{g^) we are considering here. Dropping terms suppressed by may not be numerically 
justified, but wc expect their effect to be small, in particular because this only concerns 
four out of more than a hundred degrees of freedom in the standard model. The practical 
reason is that we want to avoid computing all three-loop diagrams in the effective theory, 
as well as the renormalization and scale dependence of the Higgs mass at two-loop level. 

Evaluating all the three-loop vacuum diagrams of theory 5'ei and two-loop diagrams of 
5e2 (there are only three of them), and setting the number of fermion families to np = 3, 
we get 



PM2 _ dp ~3 _ 
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, (3.20) 
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The 0{g^) two-loop divergences in pMi cancel against as before. In addition, there 
are 1/e terms left in pE whose coefficients by themselves are of order (7^ but combine to 
a term proportional to m|, which cancels against the divergence in pm2- Unlike in the 
high-temperature case, there are also 0{g^) divergences in pmi, but these go away when 
the one-loop corrections in m| multiplying the "sunset" diagram in pM2 are included. The 
scale dependence cancels along with divergences, if we take into account the running of the 
couplings in the full theory, and set all regularization scales to be the same, A = //3 = /is. 
The phase transition can now be safely approached, since the result (I3.20p for pM2 is 
perfectly well-behaved as mi goes to zero. 



3.3 Numerical results 

To see how the multitude of terms we have computed affects the physical pressure, we 
have to supply some numbers for the parameters of the 4d Lagrangian, Eq. (|3.ip . Apart 
from the yet undiscovered Higgs particle mass, the standard model parameters have been 
measured to great precision in collider experiments, in particular LEP. The values of 
couplings can be determined from their tree-level relations to various mass parameters, 

g'^imz) = 4\/2G^m^ , 
g^imz) = 2V2G^ml, 

where m\Y = 80.40 GeV, raz = 91.19 GeV and mt = 174 GeV are the masses of the W 
and Z bosons and the top quark, respectively, and = 1.664 • 10~^ GeV~^ is the Fermi 
coupling constant [60]. The cited values are what we have used in [2,3], and they remain 
practically unchanged in the more recent Review of Particle Physics [61], with only the 
strong coupling being slightly smaller, as{mz) = 0.1176. Searches for the Higgs particle 
give its mass a lower limit mn ^114 GeV but leave it otherwise unknown. We have used 
the value mn = 130 GeV in all our analysis. The pressure of the full standard model is 
very insensitive to the Higgs mass when we are not close to the phase transition, although 
Tc itself depends on mn- In [2] we have shown that increasing rriH to 200 GeV causes a 
relative change of 10"'^ in the pressure. 

In addition to the standard model, we have also studied the simpler SU(2) + Higgs 
theory, for which the corresponding results can be found from those computed above by 
setting g'^ = g^ = gy = np = 0- Besides simpler analytic expressions, this model has some 
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Figure 3.2: Pressure of the standard model at different orders of perturbation theory. 

advantages over the standard model when we want to study the behavior of dimensional 
reduction near the critical temperature. The Higgs field drives the phase transition, but 
it only represents four of the 106.75 effective degrees of freedom in the standard model, 
so its effects on the pressure are hard to see. In SU(2) + Higgs model the corresponding 
number is only 10. Also, this toy model is weakly coupled, with A ~ 0.20 and ~ 0.61 
given by Eq. (|3.22|) . so the perturbative expansion converges better than in the presence 
of large couplings Qs and gy- More details on computations in this model can be found 
in [62]. Thermodynamics of SU(2) + fundamental Higgs theory have also been studied on 
lattice [63]. 

In Fig. 13.21 we have plotted the pressure of the full standard model up to 0{g^), com- 
bining the results in Eqs. ()3.7p . ()3.16p and the pure QCD pressure taken from [9-12] (with 
one- loop quark diagrams subtracted to avoid double counting), 

p{T) = pe{T) + pm{T) + pqcd(T) + ■ 0{g^). (3.23) 

The radiative electroweak corrections also need to be taken into account in the two-loop 
electric gluon mass m|g which we insert in Pqcd(^) above. They have been computed 
in [2], and used in the leading order term 2m^^/'i-K. 

Although the physical pressure is independent of the renormalization scale, the per- 
turbative expansion depends on the scale through the renormalization of parameters at 
orders higher than those included in the computation, and we need to fix the scale to define 
the couplings. As the remaining scale dependence is cancelled by higher order terms, the 
magnitude of the unknown corrections can be estimated by varying the scale. We have 
chosen to use A = 27rT, having shown that the scale dependence is indeed weak. 

The pressure in all our plots is normalized to the Stefan-Boltzmann result of non- 



25 



interacting gas of relativistic particles, 
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106.75 
10 



(3.24) 



for standard model and SU(2) + Higgs, respectively. The number multiplying T*^ is 
actually uei of Eq. (j3.7p + the gluon contribution. 

Fig. 13.21 shows that the perturbative expansion does not converge very well at moderate 
temperatures, but instead the 0{g^) correction is even larger than any of the preceeding 
terms. This behavior is known in QCD, and the strong coupling constant is still large 
at electroweak temperatures, gl{mz) ~ 1.48. In the full standard model the strongly 
interacting degrees of freedom sum up to 79, or 74% of the number in Eq. (j3.24p . so the 
terms with gluon exchange clearly dominate the higher order corrections, together with 
the few terms containing Yukawa interactions. The pressure lies 5-10% below the ideal 
gas result, and begins to converge very slowly at temperatures in TeV range. 

Close to the phase transition we have to use the resummed result of Eq. (j3.17p for the 
pressure. In Fig. 13.31 we have plotted both the high-temperature result pwt{T) and light 
scalar mass resummed result p-p"£{T). As the figure shows, the high-temperature result is 
very sharply peaked at Tc, whereas the corrected computation goes through the transition 
smoothly. Of course, below Tc the system goes to the nonsymmetric ground state whose 
pressure is larger than the symmetric phase pressure plotted in Fig. 13.31 since a thermal 
system always tries to minimize the free energy, or maximize the pressure. When m| 
becomes negative slightly below Tc, the symmetric phase pressure develops an imaginary 
part, which can be interpreted as the decay rate of the unstable symmetric phase [64]. In 
this region we have plotted the real part of the pressure in our figures. 

The two curves in Fig. 13. 31 differ by a term roughly proportional to T^, but the difference 
is only about 0.3%. This is because in ppx we have resummed another class of diagrams 
that are not present in pwr- We have also left out all three- loop diagrams in S'e2, which 
would contribute at order rh'^g^, and the 0{g^) corrections to mg. In particular the terms 
with strong and Yukawa couplings might be important even at this order. 

Corresponding plots for the simpler SU(2) + fundamental Higgs theory are shown in 
Fig. 13.41 and Fig. 13.51 with parameters taken from Eq. ()3.22p using mw = 80 GeV and 
niH = 130 GeV. In both figures we have also plotted an approximation of the broken phase 
pressure below Tc to indicate the phase transition at Tc ~ 220 GeV. It can be derived from 
the two- loop (order g^) computations of the effective potential [54] and the pressure in 
the symmetric phase by using pBp{T,(l)) = psp{T) — Ves{T,(f)). Fig. 13.41 should not be 
trusted near the transition, but away from Tc it shows that in the weakly coupled model 
the perturbative expansion converges nicely and settles to a level which is about 2% below 
the ideal gas pressure. Note that in the one-loop renormalization A(27rT) grows with T, 
so there is no reason to expect free theory behavior even at very high temperatures. 

Close to the phase transition Fig. 13.51 shows a similar peak in pht as in the standard 
model, while ppt does not see the transition at all until mg becomes negative at about 
15 GeV below Tc- The constant (times T^) difference between the two pressures is larger 
than in the standard model, about 2%. The couplings are all small here, so we expect that 
the three- loop diagrams in in 5e2 and the 0{g'^) corrections to the fundamental scalar mass 
are not important. However, in this model larger fraction of the overall pressure comes 
from the Higgs sector, so the result is more sensitive to different resummations on the 
scalar propagator. 
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Figure 3.4: The pressure of SU(2) + Higgs theory. 
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Figure 3.5: The pressure of SU(2) + Higgs theory near the phase transition. Here pht stands for 
the high-tcmpcraturc calculation, while ppx lias Higgs mass corrections resummed. Also shown is 
the 0{g^) pressure in the broken symmetry phase. 



We have computed the pressure of the full standard model to order both near the 
phase transition and at high temperatures. In principle, it is possible to go one step further 
and compute the coefficient of the last perturbatively accessible term of order In g. This 
computation has been carried out in QCD [12], but because the QCD pressure dominates 
the standard model pressure and gg is large, we expect advances in understanding the QCD 
pressure to be more important than four-loop diagrams in electroweak theory. Another 
computable term would be the neglected three- loop diagrams and 0{g*) mass corrections 
in S'e2 close to the phase transition, contributing parametrically at order g^'^ ^ but these 
affect only the Higgs sector and are probably too small to have any physical implications. 
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Chapter 4 

Mesonic correlation lengths 



Dimensional reduction is only useful in cases where the physical quantity of interest is time- 
independent. We can then average over any time coordinates, or, in terms of momentum 
space Green's functions, take the limit pQ = uj ^ on all external momenta. This is 
trivially so in the case of thermodynamic potentials computed in the previous chapter, 
since they are computed from vacuum diagrams, having no external legs at all. 

Another class of observables undergoing a dimensional reduction are various screening 
correlators, which describe the response of the system to a time-independent external 
perturbation. At low momenta they are usually dominated by simple imaginary poles in 
the momentum space, leading to spatial correlators that at large scales decay exponentially 
with the distance. The characteristic scale of this exponential fall-off is referred to as the 
screening length and its inverse the screening mass. A typical example would be the electric 
field of a point charge immersed in electromagnetic plasma, which at large distances is 
screened by the electric mass m^j ~ It should be noted that in general the masses 

of real-time bound states can be very different from the corresponding screening masses. 

Static correlators of bosonic operators have been succesfully studied using dimensionally 
reduced effective theories. For example, various gluonic correlation lengths have been 
measured by implementing the three-dimensional theory of static gluons (EQCD) on lattice 
[15-20]. In these works the fermionic modes are integrated out as in sections 12.21 and 13.11 
leaving a theory of soft gluonic excitations around the perturbative vacuum. However, this 
is not the only option when deriving an effective theory, but we can as well expand around 
any other saddle point of the action, corresponding to a choice 7^ in Eq. p. 170 . In 
particular, not all fermionic modes need to be integrated out, but the effective theory can 
constructed around some specific fermionic state. This opens a possibility to study also 
fermionic correlators using dimensional reduction. 

Of particular interest are operators consisting of a light quark-antiquark pair propa- 
gating in the hot medium. In [65] it was suggested that at scales comparable to the 
magnetic scale l/g^T the spectrum of quark-gluon plasma consists of color-singlet modes 
only, while the colored excitations are dynamically confined. The lowest lying excitations 
at these scales would then be the various glueball modes and the mesonic and baryonic 
states consisting of two and three quarks, respectively. In order to better understand the 
long-distance behavior of the plasma, the properties of these states have been measured 
in detail on lattice. Most of these studies have been devoted to Euclidean correlators, 
or the screening properties of these operators, due to the inherently Euclidean nature of 
lattice simulations. In particular the spectrum of the hadronic screening masses has been 
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carefully measured, the first simulations dating back 20 years [66,67]. 

While combining perturbative calculations with lattice simulations using dimensional 
reduction has been very useful when measuring the glueball spectrum, the hadronic screen- 
ing masses have been measured using expensive 4d simulations. However, the need for 
analytical tools is even greater in the fermionic sector, where the lattice simulations have 
difficulties in treating the light dynamical quarks correctly. The situation is yet worse 
when we allow for quark chemical potentials, which make the fermion determinant com- 
plex and ruin the conventional importance sampling. On the other hand, operators built 
out of quark fields are usually less infrared sensitive, so perturbation theory should more 
applicable in computing their properties. 

The first attempts to determine the screening masses of mesonic states at high tempera- 
tures using dimensional reduction were more of a qualitative nature, since they knowingly 
left out corrections of the same order as the leading term [68,69]. In particular, the scale 
inside the logarithm in the two-dimensional Coulomb potential ~ Inr was not fully iden- 
tified. A more systematic approach was developed by Huang and Lissia in [30], where it 
was shown that the dimensionally reduced theory for fermionic modes can be formulated 
in terms of massive non-relativistic quarks in 2-1-1 dimensions. They also discussed the 
correct power counting of different operators, and computed one-loop corrections to the 
quark self-energy and the quark-gluon interaction vertex. However, although the efi^ective 
theory was derived in order to calculate screening quantities, the authors did not proceed 
to compute any masses in that work. 

Following [30], wc used similar methods in [1] to derive an cff'cctive three-dimensional 
theory for the lowest fermionic modes zLttT, which dominate the mesonic correlator at 
large distances. This theory takes the form of non-relativistic quarks coupled to EQCD, 
or the spatial gluons and an adjoint scalar field. Because the fermionic sector of the 
theory is very similar in form and power counting to the effective theory for heavy quarks 
in four dimensions known as "non-relativistic QCD", we have named the reduced theory 
NRQCD3. Using this theory, we were able to compute the next-to-leading order correction 
to mesonic screening masses in perturbation theory, and this computation was extended 
to finite quark chemical potentials in [4]. In this chapter we will review these results. 

4.1 Linear response theory and screening phenomena 

The correlation functions usually computed in theoretical calculations are related to phys- 
ically measurable quantities through linear response theory. Our short presentation here 
follows mostly [36], and is somewhat biased towards screening physics. 

Consider perturbing the system in equilibrium with some external probe, described by 
an interaction Hamiltonian V{t) which vanishes for t < 0. In the Schrodinger picture the 
time-development of an unperturbed state is given by the time-independent Hamiltonian 
H, while the effect of V{t) can be written in terms of a time-development operator U{t), 




iVs(O)) = e 

umsm 




= H\i;'s{t)) + e-'''HdtU{tMsm = (H + VM'sit)) 

idtU{t) = e'"*V{t)e-'^*U{t), U{t) = 1 for t<0. 



(4.1) 
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In the last equation Vnit), the potential in the unperturbed Heisenberg picture is seen 
to emerge. If V is small, U{t) can be solved recursively as a series in V by integrating 
Eq. (BH), 

Uit) = 1 - Tdti VHih) - fdh Tdta VHiti)VH(,t2) + OiV^). (4.2) 
^0 Jo Jo 

The change in the expectation value of an arbitrary operator 0{t) in the Schrodinger 
picture is then 

5{d{t)) ^ {^'s{t)\d{tWs{t)) - mt)\d{t)\Mt)) 

POO 

= -i At'{M0{t-t')[dH{t)Mt')]\'4^H). (4.3) 



JO 

where the operators and the state vectors are now all in the Heisenberg picture with the 
unperturbed Hamiltonian H. In particular, Eq. (j4.3p applies to the eigenstates of H, 
so we can sum over all states in the ensemble with appropriate weights, and replace the 
expectation value in a specific state by thermal average. 

Typically, the external interaction can be written as a time-dependent c-number source 
v{t, x) coupled to the system through some current J built of field operators, 

V{t) = jd?xJm,^))v{t,^), (4.4) 

and the response of the system is measured through the same current. Eq. (|4.3|) can then 
be written in terms of the retarded correlation function Dr^ 



5{j{4>{t,^))) = -i I dt' j<5?x'v{t',^){e{t-t')[j{<t>H{t,^)).J{4>H{t',^))]) 

dVL>/j(x -x')w(x'), (4.5) 



where the lower limit of the time integration can be extended to — oo because f (t, x) = 
for t < 0. In this expression the response of the system is clearly separated into the 
retarded propagator, which is specific to the thermal system in question, convoluted with 
a factor v{x) depending on the details of the perturbation. 

The excitations of the system manifest themselves as large responses to an external 
perturbation. Going into the momentum space, the Fourier transform of Eq. ()4.5p reads 



5{J){uj, k) = -iDniuj, \l)v{uj, k), (4.6) 

so this is equivalent to having a pole in the propagator for some values of frequency uj and 
momentum k. In general there are some real-time excitations, whose frequencies depend 
on the momentum through the dispersion relation uo = cj(k). In addition, at the static 
limit a; ^ the propagator may contain poles at imaginary momenta, corresponding to 
perturbations that decay exponentially with distance due to plasma screening. The former 
are difficult to study in imaginary time formalism because of the need to analytically 
continue the imaginary time results, but the time-independent screening correlators are 
well-adapted for Euclidean computations. 

It should be noted that the correlators often have very different dependence on oj and 
k, since the thermal environment breaks the Lorentz invariance. The screening states 
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and the real-time excitations describe completely different physics, and the corresponding 
poles need not be related at finite temperature. One should also be careful when taking 
limits k — > 0, since different orders of limits may give different results, and the correct 
procedure depends on the physical situation. A good comparison of screening and real- 
time quantities in the context of solvable 2+1-dimensional Gross-Neveu model can be 
found in [70]. 

4.2 Mesonic screening states at high temperatures 

The Euclidean Lagrangian governing the behavior of quarks and gluons at finite temper- 
ature is 

£e = + + M)i; , (4.7) 

where the covariant derivative is defined as D^ip = d/^iip — igA'^T°"il: and the gluon field 
strength Ffj_y = i/g[D^,D^], being the generators in the fundamental representation 
of SU(A'^c)- We consider Np flavors of degenerate quarks, so the quark field ip is an Np- 
component vector in flavor space, and the mass matrix M is proportional to unit matrix, 
M = m ■ Ij^p. For simplicity we will set m = in most of what follows. 

Finite quark densities are included through chemical potentials /if multiplying the quark 
number density operators A/f = tpfjoipi- This is precisely the same structure as in the 
momentum time-component, so the actual effect when doing computations in perturbation 
theory is a shift in pQ, which in absence of chemical potential woud be one of the fermionic 
Matsubara frequencies, i[{2n + l)7rT] i[{2n + 1)ttT — i/if]. In pure QCD all quark 
numbers are conserved separately, so the chemical potential for each flavor can be chosen 
independently. Weak interactions, on the other hand, mix different quark flavors, mostly 
inside SU(2) doublets but also between families, because the weak interaction eigenstates 
differ from the mass eigenstates (for a brief review on this mixing, see [61]). Because of 
this, only certain combinations of the baryon number and the different lepton numbers 
are conserved in the full standard model. We assume that the time scales of chemical 
equilibration through weak interactions are much larger than the characteristic time scales 
of QCD processes we are studying, even when discussing static correlators, and continue 
to use independent chemical potentials for ech flavor. In the numerical studies we will use 
two distinct cases, isoscalar (/^u = /^d = fJ-s) and isovector (/Xu = —fi^ = ^y) chemical 
potentials for illustration, but the analytical results are applicable to general /if. 

The quark fields il),ip can be used to define mesonic operators of different spin and flavor 
structures. These operators can be thought of either as the currents coupling to external 
perturbations as in Eq. ()4.4p or as interpolating operators for physical particle states, 
although in the latter case it should be remembered that the real-time mesonic bound 
states do not survive at very high temperatures, and that their properties may be very 
different from the corresponding screening states. We denote 

= ^jF^Tip, (4.8) 

where P is one of {1, 75, 7/^, 7/^75} for scalar, pseudoscalar, vector and axial vector objects 
Qa _ S"',P"',V^,A'^, respectively. The flavor strucure is written in terms of the identity 
matrix F'^ and the traceless matrices F"", which satisfy 

F' = 1n^, TT[F''F'] = h''\ a,6 = l,...,iV2 -1. (4.9) 
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(a) (b) (c) (d) 

Figure 4.1: Classes of diagrams contributing to the meson correlator. Diagrams of type (d) only 
apply to flavor singlets. 



We wish to compute the static correlators of the above operators. The particular cor- 
relators we are interested in are defined as 

fl/T f 

C^[0'',0^] = j dT J d^xe^'^-''{O''{T,x)O''{0,0)). (4.10) 

The quantity in this equation seems different from the retarded propagator defined in 
Eq. (j4.5p . The number of relevant propagator-like functions at finite temperature is large 
because of the possibilities of having either real or imaginary time coordinates as well 
as different orderings of the operators, but they are all related to the spectral function 
p{io,q), which in turn is given by the analytic continuation of the Euclidean correlator in 
Eq. ()4.10p . All these relations can be found in standard textbooks (see e.g. [36]), and a 
good summary is given in [71]. The relevant relation for our computations is 

/ ON f°° dujp{uj,q) /■°° duj p{uj,q) ,. ^ / x 

DR{qo = 0,q)= — = lini / — ' ^ = lim Cs(go,q), (4.11) 



where CEiqo^^) is the Euclidean correlator in Eq. ()4.10p with an arbitrary (Euclidean) 
momentum zero-component go- 

We can use the rotational invariance of the system to choose the vector x, the direction 
in which we are measuring the correlations, to point in the direction, and further 
average over the transverse xij;2-plane. The fundamental quantity we are studying here 
is then the z-dependent correlator 

C,[0^0^] = J dVC(xx,.)[0",0'] = J^^^dT J dV(O'^(r,xx,z)O^0,0,0)). (4.12) 

At very high temperature the QCD coupling is small due to asymptotic freedom, and 
the correlator Cz can be computed using perturbation theory. The leading order result 
is given by the free theory diagram Fig. 14.1( a) consisting of two noninteracting quarks 
propagating in the hot medium. When the chemical potentials all vanish, this diagram 
can be computed in the momentum space, and apart from constants corresponding to 
terms ~ 6{x) in the coordinate space, the result is proportional to the function 

0/0 N_ fdh> 1 i 1 - iq 

^3d(2-„) = J (^[,2+p2][,2 + (p + q)2] - ^1-^Z;;:T^ ' (4.13) 

summed over all fermionic momenta with coefficients that depend on the spin structure 
of the operator in question. At large distances, or low momenta, the behavior of the 
correlator is dominated by the lowest singularities of this expression, located at q = zizi2irT. 
If we for a moment imagine rotating the three-dimensional theory of the lowest Matsubara 
modes to 2-1-1 dimensions and treating the direction of q, or X3, as a time coordinate, 
this singularity becomes a branch cut on the real ga-axis, and can be understood as the 
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threshold of producing two free quarks of mass vrT. When the interactions are turned on, 
we expect this singularity to convert into a pole corresponding to a bound state in the 
2+1-dimensional theory. 

It should be noted that the temporal components of vector and axial vector correlators 
do not have the singularity associated with B^^{2uJn), since carrying out the Dirac algebra 
gives these terms a prefactor + 4u;^, which regularizes the singularity at q = ±2iuJn- 
Computing in the configuration space, it is easy to see that these correlators are suppressed 
by powers of distance and decay even faster, since the contribution from the lowest sin- 
gularity is removed. The longitudinal components V3 and A3 of those correlators vanish 
completely, apart from contact terms, because of the current conservation 9^V^ = 0. 

When the chemical potentials are turned on, even the free correlator becomes very hard 
to compute. This is due to the shift in the temporal momentum components, which causes 
the correlator to mix different Matsubara modes even after integrating over r-direction. 
If the two quarks have identical chemical potentials, which is the case for flavor singlets 
or isoscalar chemical potential, the correlator Cz in Eq. (I4.12|) can be computed explicitly. 
For scalar operator the result is 



NrT 



Sttz sinh 2itTz 



2ttT coth 2'kTz cos 2^z + 2^ sin 2/iz H — cos 2^z ) (4.14) 



87r2z3 



1 



+ — '^ + —-^]{'^T^Tzf + Oiz^^ 
V360 6 7r2r2 12 7r4r4 ' ^ ' ^ 



^ab ^cT'^ 1 c~^"^^'^ 



1 \ /i ■ 

1 + ^^j^ I COS 2^iz + — sm 2nz 



+ Oie 



-iTTTz\ 



where we have also indicated the limiting behavior at small and large distances, respec- 
tively. For ;U = this result agrees with the previous zero density computation in [72] . The 
effect of a small chemical potential on the correlator is seen as oscillations with wavelength 

= K I [i inside the zero density envelope, while for a large [i S> vrT the interaction with 
the particle bath is so strong that correlator oscillates wildly, averaging to zero. 

The case of arbitrary chemical potentials is much harder because of the summations 
mixing different modes, but a good approximation can be found by only taking into account 
the lowest modes = ibvrT, which dominate the correlator at large distances. The scalar 
correlator can then be written in terms of exponential integral function Ei(z), which we 
approximate at large z to leading order in l/z, resulting in 



/V 7^2 



27r 
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A^ij (27rr)2 + Ji 



2 ^ij sm/XijZ 



+ 



A/iij 



+ 



A^i: 



(2vrr)2 + /i2^ 



27rrcos /iijZ 



(4.15) 



where A/iy = /i; — /ij and /iij = /^i + /.ij. Details of the computation and a more complete 
result for the general case can be found in the Appendix A of [4]. In this correlator 
the scale of oscillations is /iij, a direct generalization of the isoscalar case. In particular, 
for isovector chemical potentials the oscillation vanishes (although it reappears at terms 
suppressed by Xjz as sinA/iyz, see [4]) and the correlator has the same functional form 
as in the = case. 

Together all these results show that the free correlator falls off as exp(— 27rr2;), regardless 
of chemical potentials. The finite density only shifts the singularity in ^3 by an imaginary 
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Figure 4.2: The free scalar meson correlator for different values of isoscalar chemical potential. 

part, which shows as oscillations in Cz- We define the screening mass as the coefficient of 
this exponential fall-off, or the real part of the pole location, because that is what deter- 
mines the asymptotic behavior of the correlator. In Fig. (14. 2p we have plotted the isoscalar 
correlator Eq. (j4.14p for different values of chemical potential. The figure shows that when 
the wavelength of the oscillations = vr/// is large, /i < O.SttT, the oscillatory behavior 
becomes apparent only at distances where the correlator is exponentially suppressed, and 
it is hard to discern the cosine term from an increased exponential fall-off. For larger ^ 
this difference is obvious, and even more so when studying the asymptotic behavior of 
the correlator, so in the leading order of perturbation theory m = 2'kT + 0{g^) is the 
consistent definition also at finite density. 

4.3 Effective theory for the lowest fermionic modes 

Going beyond the leading order, the location of the singularity in the static meson corre- 
lator is modified by interactions with gluons. We would like to compute in perturbation 
theory the next-to-leading order corrections to the screening masses, which requires com- 
puting the diagrams of types Fig. I4.1l fb.c). Note that the flavor singlet correlations are 
also mediated by purely gluonic states, some of which have masses lower than 2'nT. These 
states are completely non-perturbative, so we only concentrate in this work on flavor non- 
singlets. The couplings of different quark operators to glueballs have been worked out in [1] 
and the masses of these gluonic operators for both zero and finite chemical potentials have 
been measured in [19]. 

Computing the one-gluon diagrams is not enough, however, since one runs here into 
the same infrared problems as always when computing with light particles at finite tem- 
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perature. In particular, for soft momenta (which imphcitly requires bosonic zero mode 
iOn = 0) the one-loop correction to the gluonic zero mode propagator is of the same order 
as the leading term and needs to be resummed into the propagator to get rid of infrared 
singularities arising from soft gluons. As discussed in previous chapters, at order this 
resummation gives the static temporal gluon component an electric mass rn-E ~ gT, which 
we have to include in the soft gluon propagators of diagrams Fig. I4.1( b.c). 

Another class of diagrams requiring resummations consists of the graphs with soft gluon 
exchanges, as shown in Fig. l4.1( b). The integration over when going to coordinate space 
gets the largest contribution from the poles at ~ ibzvrT, and if the gluon momentum is 
small, the additional quark and gluon propagators in the diagram are also nearly on-shell, 
1/^ ~ 0{1/ g'^T), compensating for the factors of g from the vertices. To get a consistent 
next-to-leading order correction we therefore have to sum over all diagrams with an arbi- 
trary number of soft gluon exchanges. This is a common requirement for having a bound 
state in theory, since summing only a finite number of diagrams with free propagators 
cannot give rise to new singularities, and it applies as well to the screening states under 
discussion. We wish to avoid the exceedingly complicated formalism of relativistic bound 
states, and instead make use of the hierarchy between the Matsubara modes vrT and the 
momentum scales where the resummations become necessary, p ~ gT. This suggests using 
an effective theory to organize the resummations of the soft modes, while simply including 
the effects of the large momenta in the parameters. 

The leading order computation shows that at large distances the screening correlator 
is dominated by the lowest Matsubara modes zbvrT, while the contribution coming from 
the other modes is exponentially suppressed. We concentrate on the lowest modes and 
integrate over all excitations of momenta ~ vrT around these modes, which gives an 
effective theory only containing the static (a;„ = 0) gluons and the fermionic modes with 
ujn = ±7rT. In particular, all the other fermionic states are off-shell by 2tiT ^ gT, so this 
theory does not contain any creation or annihilation of quarks, but the only quark lines are 
those entering and leaving the diagrams as external legs. Because we are computing near 
the screening pole, these external quarks are almost on-shell, = 0, and the interactions 
with soft gluons do not change that very much. The relevant expansion parameter is 
then the "off-shellness" Pq -|- +P3 ~ Q^T'^-, which is of the same order as the momenta 
of the soft gluons. If we rotate the 3-dimensional action of the quark modes into the 
2-|-l-dimensional Minkowski space, the Matsubara mode ztTrT can be viewed as a heavy 
quark mass, and the restriction to momenta much lower than this and the separation of 
quarks from antiquarks effectively makes the quarks non-relativistic in 2+1 dimensions. 
Note that there is nothing special in the lowest Matsubara modes, but we could equally 
well derive a similar theory for modes w„ = SvrT, for example. The reason we have singled 
out the lowest modes is because they dominate the screening correlator, which we aim to 
compute. 

The difference between the purely gluonic EQCD and the effective theory considered 
here is the state around which the expansion takes place, and is dictated by the physical 
application we have in mind. EQCD is directed at computing either vacuum diagrams 
or gluonic correlators, which do not have any quarks on external legs. Any fermionic 
excitation is then off-shell at least by 27rT, and only contributes at the high-momentum 
integration. The mesonic operators, on the other hand, are built out of quark fields, which 
cannot be integrated out if we intend to compute with them. Still, given an external quark 
state, any additional quarks would again be very much off-shell, so expanding around that 
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state in low momenta we can neglect all other fermionic excitations. To compare with 
QED, the relation between the two theories roughly corresponds to the Euler-Heisenberg 
effective Lagrangian for photon-photon interactions and the non-relativistic quantum me- 
chanics used to compute hydrogen (or positronium) binding energies. In the former the 
electrons have been integrated out completely, whereas in the latter we can mostly ignore 
the contribution of the states with e.g. two electrons and a positron, not because the mass 
would be much higher (3me vs. me), but because they are off-shell by ~ 2me, which is 
large compared to typical momenta ~ airie- 

The bosonic sector of the dimensionally reduced QCD is well-known [5], and in finite 
density it reads [19] 



= ^Tr 4 + Tr [A, Aof + m|Tr^2 + ^ ^^^^^3 + (Tr ^g)^ ^ ^(2)^^^! , 

f 

(4.16) 

where we have also included the cubic Aq self-coupling proportional to quark chemical 
potentials. The parameters are found by matching gluonic 2- and 3-point functions. We 
only need the adjoint scalar mass to one-loop level and other couplings at tree level, so 

m| = 5^T2(^^ + ^ + ^E^)' 9l=9'T, A(^'2) = 0(5^T). (4.17) 

The cubic and quartic Aq self-interactions can be ignored, since they contribute to meson 
correlators only at order and higher. The bosonic theory in Eq. (14.160 describes a 
3-dimensional gauge theory with a massive adjoint scalar ^q, and has been extensively 
used to compute gluonic quantities at high temperatures, as discussed in the beginning of 
this chapter. 

On the fermionic sector the tree- level Lagrangian is just a sum over the parts of the full 
theory Lagrangian containing modes a;„ = ibvrT. For a single mode w„ this term reads 

Cq = i> [i-fouJn + 7oAt - ^570^0 + IkDk + 73-D3] ■0, (4.18) 

where k = 1,2 and Aq is the gluonic zero mode. We have separated the xs-direction 
from the other spatial components, anticipating the choice to measure correlations in that 
direction. Note that the interaction with static gluons does not mix different fermion 
Matsubara modes, so we have a separate term like Eq. (14.180 for each mode we wish to 
compute with. 

Using a non-standard representation for Dirac matrices (see [1,4] for details) and de- 
composing the four-component spinor as 

^P=(l], (4.19) 



the Lagrangian can be written in a form where the fields x a-^d (p are light and heavy close 
to the pole = ipo = i{oJn — «^), respectively, while the roles are reversed at the other 
pole p3 = —ipo- Solving the equation of motion for the heavy component and expanding 
the resulting non-local operators in 1/poi we get the non-relativistic Lagrangian 



£q ^ ix^ po - gAo + ^3 - ^ (^1 + ^.[(^i, 

po - 5^0 - -D3 - ^ [dI + (rj]F, 



X 



+ 0{-^J. (4.20) 
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Figure 4.3: One-loop correction to the quark self-energy. 



All dependence on the chemical potential at this level is contained in the shift of the 
temporal momentum component, pQ = cOn — ifJ-- The Lagrangian is easier to understand if 
we again imagine rotating to 2-1-1 dimensions and setting z = it. The zero-point energy 
is then given by po, while the other free terms combine to —{idt + V^/2po)j the standard 
nonrelativistic kinetic term with mass pq. If we forget the relativistic origin of these terms, 
there is no reason for the zero-point energy and the mass parameter in the kinetic term 
to be the same, and as we will see shortly, the loop corrections will give these parameters 
different values. It should be noted that for a;„ > the field (p has negative mass, and 
should be interpreted as the antiparticle of x- 

Already at the tree-level the expansion in 1 /pQ gives rise to an infinite number of terms, 
and beyond this we will have to take into account all possible terms allowed by symme- 
tries. To limit the possibilities, a power counting has to be established. Requiring all the 
transverse momenta to be at most of the order of the electric mass, p_L < gT, and that 
the terms in the action to be of order unity, we get 



X 



l/|xx| ~ (x3/xi)i/2 ^ ^1/2^1/2 _ (4_21) 



The off-shellness Ap3 = p^, Upo corresponds to the kinetic energy in the 2-1-1 dimensional 
theory, as can be verified from the poles in the quark propagators, Eqs. (I4.25p . ()4.26p . For 
nearly on-shell quarks we can then estimate the derivative by 



Ap3 ~ p±/po ^ 9 T =^ ds g T acting on quarks. (4.22) 

Using this power counting, and keeping only terms required to give the screening mass to 
order g'^, the final form of the fermionic Lagrangian is 

4ff = iX^ [m - gEAo + D3- ^) X + - se^o - ^3 - ^ , (4.23) 

where only the zero-point energy, which we will henceforth denote by M, needs to be 
matched beyond tree-level. This matching is carried out by comparing the poles in the 
one-loop corrected quark propagator depicted in Fig. 14.31 The details of the computation 
are given in the original papers, and the result is simply 

M = po + g^Cp^ (l + =ujn-if^ + g^Cp—^^— (1 + . (4.24) 



8po \ 7r2r2 ; " ' y 8(^^ _ -^^ ^22-2 J ■ 

It should be noted that for Un = ±vrT, which are the modes of interest here, the real part 
of M does not depend on fj.. 

The free quark propagators following from Eq. (|4.23p are 

{Xuip)xl{q)) = Suv{2nf6{p - q) ~' (4.25) 

M + ip^ + Y>J2po 

{MpWv{q)) = Suv{2nf5{p - q) ■ 2 /n (4-26) 

M -ip3 + py 2po 
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or in the configuration space 



{Xu{x)xl{y)) = -iSuvd{uJn{x3 

{(pu{x)(l)l{y)) = -iSuvO{u;n{y3 

In these equations it is obvious that for a;^ > the field x propagates forward and 
(p backward in X3, the time coordinate of the 2+1-dimensional theory, confirming our 
interpretation of 4> as the antiparticle of x- For negative modes these roles are reversed, 
with (j) propagating forward in X3. In [1] we expanded these propagators in p_l/po inside 
loop integrals to make sure that the transverse momenta are parametrically smaller than 
the heavy scale ~ T. However, as discussed in [73], while this works well for (single) 
heavy quark effective theory, in NRQCD there are difficulties at two-loop level with this 
approach, and it is preferable to keep the kinetic terms summed into the propagators. 
For our modest purposes there is no real difference, but in [4] we chose not expand in 
transverse momenta. For consistency, this requires that we expand the gluon field in 
multipole expansion, which at this level boils down to disallowing transverse momentum 
transfer from gluons to quarks. The masses and the quark- antiquark potential turn out 
to not depend on the way we treat the propagators, whereas in order to compute the full 
correlators it is necessary to keep the kinetic terms resummed. As a check, we have shown 
that NRQCD3 is able to reproduce the leading order scalar correlator in Eq. ()4.14p exactly. 

In terms of the new fields x 4' the operators whose correlators we intend to compute 
can be written as 
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(4.29) 



It should be noted that the temporal and longitudinal components of vector and axial 
vector currents consist of terms like x^X so their correlators are proportional to 

6{z)0{—z) and vanish at nonzero distances. For V3 and ^3 this follows from current conser- 
vation, whereas the correlators for charges Vq and ^0 are power-suppressed, as discussed 
in the previous section. Apart from those operators, the correlator in the effective theory 
is independent of the spin structure, up to a multiplicative constant, since (o",)^. The 
flavor structure on the other hand is significant if we allow for finite chemical potentials. 

4.4 Solving the screening states 

Having derived the Lagrangian for NRQCD3, which is just the sum of the bosonic part in 
Eq. (I4.16P and the fermionic part in Eq. ()4.23p . it remains to find the masses of mesonic 
operators. These correspond to bound states in the 2-|-l-dimensional theory, following 
from the summation of diagrams of type Fig. I4.ir b.c) with an arbitrary number of soft 
gluon exchanges. In a nonrelativistic theory the resummation can be carried out by finding 
the static potential for a x*4' pair, and then solving the resulting Schrodinger equation 
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(4.27) 
(4.28) 



39 



with this potential. Here "static" should be understood from the 2+1-dimensional point 
of view, which means that the potential will be valid for large X3. All these computations 
can be performed using the effective theory just derived, since it is only the iVn = gluons 
with low momenta that need to be summed beyond the leading order. 
The static potential is written as an expansion in g'^r, 

V{r)r^ gllnr + g^r + 0{g^r^). (4.30) 

Using either the power counting 1/r ~ p_L < gT or the leading order Schrodinger equa- 
tion we see that g'^r ~ g. The leading Coulomb-type term g'^ In r in the potential is 
then sufficient for computing 0{g'^) corrections to screening masses, and can be evaluated 
perturbatively by computing all one-gluon diagrams in the effective theory. The leading 
logarithmic term already gives a confining potential, so there is no qualitative difference 
in dropping the linear term, which is parametrically of order g^ . 

The static potential of a (f)*x pair is computed by inserting a point-splitting in the 
correlator to give the quarks a small spatial separation in the transverse direction, and 
finding the Schrodinger- type equation satisfied by this correlator at z — > 00 limit. The 
details of this computation can be found Appendix B of [4], and the result is 

^(r) = '^{^^^+^E- K,{m^r)) , (4.31) 

where Kq is a modified Bessel function. The result is both ultraviolet and infrared finite 
once we have resummed the gluon self-energy corrections to an electric mass mE, while 
letting TTi-E — > we would again find the infrared divergences of the original theory. The 
screening masses at this order are not sensitive to the magnetic screening of spatial gluons, 
as this would manifest itself as divergences in the potential. The potential in Eq. ()4.3ip 
depends on the quark chemical potentials only through m-E, as the leading order is only 
sensitive to the propagation of gluons in the hot medium while the quarks simply act as 
static color charges. 

The Schrodinger equation satisfied by the correlator was already found as an interme- 
diate result when computing the potential, and it reads 



±(Mi + Mj)--^V2 + y(r) 



^0 = mMi*o, (4.32) 



where we have separated the variables as 

C(r, z) = ^o(r)e-™f""^ , (4.33) 

and the it signs apply for u;„ = ztTrT, respectively. The fiavors of the quarks forming 
the meson are labelled with indices and in general the fiavor symmetry between 
different mesons is broken by the different chemical potentials. The parameters M and pq 
are generally complex when computing with finite chemical potentials, but for opposite 
modes they are related by M_ = —M^ and po- = ~Po+' ^° screening masses satisfy 
by rrifuii^- = rrif^jj^. Thus we only need to compute the masses for ujn = vrT, and 
in addition these relations guarantee that the full correlator, which is the sum over all 
separate Matsubara modes, behaves as 

a[0^0''] a J^i^^Fj^2cos[Im(mMi,ij)z - Oy] exp[-Re(mMi,ij)z], (4.34) 
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where ay is the overaU phase of the (piXj correlator. This is of the same form as the 
leading order term, with the real part of the mass parameter giving an exponential decay 
while the imaginary part contributes to a cosine-like oscillation term. In particular, this 
correlator is real-valued even though the term coming from any single mode is in general 
complex. 

The Schrodinger equation with potential Eq. ()4.3ip cannot be solved analytically, so 
we have to find the eigenvalues mf^n using numerical computations. For this, we cast the 
equation into a dimensionless form, which depends on the values of the physical parameters 
only through the dimensionless combinations p and Eq defined as 

p^i^f ^ + \ g29LEo^mf^n-M;-M^. (4.35) 

One should note in p the appearance of the reduced mass, typical of two-body problems. 
The physical screening mass is given by Eq which we solve numerically, 

Re(mMi) = Re(Mi + Mj) + gl^Re{Eo) = 2ttT + g^T^ Q + Re(^o)) , (4.36) 

where in the last equality we have used the fact that for the lowest modes the real part of 
M is independent of p. 

The numerical solution is found by assuming that the ground state is cylindrically 
symmetric, solving the behavior of ^o{r) around the origin and then integrating out to 
larger r and requiring square integrability. In zero density this is easy since both the 
wave function and Eq are real, and trying different values of Eq on the real axis with the 
condition that ^oif) vanishes at large distances is very fast. We find, with somewhat 
superfluous precision, 

( 0.16368014, p = 2/3, {Np = 0) 
Eq = < 0.38237416, p = 1/2, (iV^ = 2) (4.37) 
[ 0.46939139, p = 4/9 , [Np = 3) 

for different numbers of dynamical quarks Np. The number of flavors only enters in mE, so 
even the case Np = makes sense if the creation and annihilation of quarks is suppressed, 
as in the so-called quenched simulations on lattice. When the chemical potentials are 
turned on, the numerical computation becomes more demanding, as all parameters and 
the wave function have complex values, and the solution has to be searched for in the 
complex plane instead of limiting to the real axis. Moreover, the dependence on p/ttT 
cannot be solved analytically, but we have to repeat the process for each chemical potential 
separately. 

In Fig. 14.41 we have plotted the eigenvalues Eq for isoscalar chemical potential with 
different numbers of dynamical quarks. We have also studied the case where Np = 3, 
but only the two flavors in the measured operator have nonzero chemical potentials. This 
should more closely correspond to the situation in heavy ion collisions, where the two 
colliding nuclei have flnite up and down quark densities but vanishing net strangeness. 
The figure shows that the poles move along quadratic curves off the real axis, and the real 
part of Eq, which is the contribution to the screening mass, grows with ps for dynamical 
quarks. In the quenched case the real part decreases and becomes negative at ps ~ 0.787rT. 
At large values of the chemical potential the oscillations due to the imaginary part of 
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Figure 4.4: The eigenvalue Eq with the lowest real part for isoscalar fis = . . . 1. 



£^0 become strong and the numerical integration is unstable, so it is hard to go beyond 
~ ttT numerically. On the other hand, in deriving the effective theory using dimensional 
reduction we assume that T is larger than any other mass scale, so our results cannot be 
trusted for /i > irT. 

The chemical potential enters the dimensionless Schrodinger equation only through the 
parameter p. For the specific cases of isoscalar and isovector chemical potential that we 
study numerically, the dependence is 

(4.38) 

Prom Fig. 14.41 we can see the twofold influence of the chemical potential. When Np = 0, 
increasing the chemical potential just shifts p into more imaginary values, while its real part 
stays constant. The following increase in the absolute value of p decreases Re(£'o) like in 
the p = case, Eq. (j4.37p . For dynamical fermions this effect is more than compensated 
by the increase in niE, which raises the potential and accordingly increases the energy 
eigenvalues. 

In Fig. 14.51 we plot the real part of Eq, which apart from scaling and an additive p- 
independent constant is the same as the screening mass, see Eq. (14.360 . The isoscalar data 
is the same as in Fig. 14.41 discussed above, while for isovector chemical potential we see that 
the mass decreases with pv for Np < 3. This is in agreement with the previous discussion 
on the relation between \p\ and Eq, and Eq. (j4.38p . where p for isovector chemical potential 
increases with p for small number of dynamical fermions. 

Numerically the correction we have computed is small, even near the phase transition 
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where the couphng constant itself is large. For Np = 2 the range of screening masses for 
different values of chemical potentials fits in the interval 

r 0.227, /x^/ttT = 1.0 
Re(mfuii) 2ttT + c/| x i 0.187, /i/vrT = 0.0 (4.39) 

[ 0.166, fiv/irT = 1.0 

The effective coupling gE is estimated to be g'^/T ^ 2.2 when T ~ 2Tc [50], giving the 
next-to-leading order corrections to the screening mass of about 6-8%. Nevertheless, as 
long as the coupling is large there is no reason to expect that the next correction would 
be smaller by a factor of the same magnitude. 

While the Schrodinger equation Eq. (I4.32p cannot be solved exactly, we were able to 
find a simple approximate dependence on the parameters while writing this introductory 
part. Realizing that the modified Bessel function Ko{mEr) in the potential interpolates 
between — ln(mE?'/2) and 0, we can try to estimate the potential by 

VW»#(C.1..^ + C.). (4.40) 

where we expect 1 < Ci < 2. For real values of the dimensionless parameter p defined in 
Eq. (I4.35P its effects can be scaled into the dimensionless variables, giving for the potential 
in Eq. (001) 

^o(p) =^o(l)-^lnp. (4.41) 

In Fig. I4.6l we have plotted Eq vs. In p for the numerical data we have computed, assuming 
that the expression in Eq. ()4.4ip can be extended for complex values of p as well, if the 
branch cut is introduced on the negative real axis. As the figure shows, for isovector p 
with real p the behavior is extremely well described by Eq. (|4.4ip . and the agreement is 
also good for the complex values of p. Fitting a line to both real and imaginary parts 
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Figure 4.6: The numerical data in Figs. 14.4 1 and 14.51 parametrized by the logarithm of the dimen- 
sionless parameter p. 

separately, we get slopes -0.78 and -0.74, respectively. To desired accuracy, the screening 
masses from our computations can then be summarized by the two-parameter fit from the 
left plot in Fig. S^l 

Re(mfuii) = 27rr + 5^2^-^^(0.39 -0.78 In I /)|), (4.42) 

27r 

where p can be read from Eq. (j4.35p . 

4.5 Comparison with other results 

The static mesonic correlators considered here have also been measured in the lattice 
simulations. In zero density there is a long tradition of measuring the screening masses, 
for recent results see [74-77]. Early measurements gave large differences between the 
masses of different spin structure operators, but in recent works the general picture is that 
already at ~ 2Tc the masses come close to the ideal gas result, with p (vector) meson 
slightly larger than vr (pseudoscalar) and both lying 5-10% below the ideal gas result. 
The correction we have computed above is of the same magnitude, but with different sign. 
For all Np the /i = results we have computed lie above the free theory result 27rT, 
approaching the ideal gas slowly as as{T) gets smaller due to asymptotic freedom. 

The systematic errors in lattice simulations are related to dynamical quarks, light quark 
masses, and the difficulties in going to the infinite volume limit. In [76] the analysis 
of the infinite volume extrapolation is carried out, resulting in slightly higher masses 
than those measured earlier, but even then the screening masses are clearly below the 
ideal gas result. It should be noted that our perturbative result is above 2ttT also for 
Np = 0, so the potential difficulties with dynamical quarks cannot completely explain the 
disagreement with lattice measurements. On the other hand, because the strong coupling 
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is large near Tc, the higher order perturbative corrections can be as large as the 0{g'^) term 
we have computed, at least at temperatures within the reach of lattice measurements. At 
asymptotically high temperatures, however, the perturbative calculation should be valid, 
so we expect the screening masses to cross above the free theory result at high enough 
temperature. 

Recently the screening masses have also been computed by evaluating the meson spectral 
function in the HTL resummation scheme and determining the fall-off of the correlator 
through the spectral function [78]. This method is very orthogonal to our computations, 
but the weak coupling limit of those results is very similar, and can be summarized in 
our terms by setting E'o = in Eq. (|4.36|) . The difference can be explained by the soft 
gluon contributions that were knowingly left out in [78]. At temperatures close to the 
phase transition the meson properties have been studied analytically in the Nambu-Jona- 
Lasinio model, which gives screening masses well below the free theory result [79]. It 
should be noted that this method is valid at very different temperature region, and does 
not necessarily contradict our high-temperature results. 

At finite density the lattice measurements of hadronic screening masses have been per- 
formed only recently [80]. Simulations at nonzero chemical potentials are difficult because 
of the complex fermion determinants, so these computations were carried out by expand- 
ing the masses as Taylor series in fj, around = and measuring the derivatives up to 
second order response. All measurements in [80] were carried out using Np = 2 flavors 
of staggered fermions. The leading order term is just the /u = mass which behaves as 
described above, and the first derivative vanishes for both isovector and isoscalar chemical 
potentials. For the latter this follows from symmetry properties, whereas the response 
with isovector fj, is explicitly measured to be zero. In our perturbative calculation the 
first derivative at /i = vanishes as well for both chemical potentials, as required by the 
symmetry Eo{-n) = ^o(m)- 

The second derivatives give the ^-dependence of the masses as measured on lattice. For 
isoscalar // the measured second order response rises steeply at Tc and settles somewhat 
below 2/T at higher temperatures, for both tt and p. In the isovector case the response is 
small and negative, and approaches zero as temperature is raised. At small these results 
are in qualitative agreement with our perturbative calculation, but for isoscalar the actual 
numbers differ by orders of magnitude. For comparison, we have fitted quadratic curves 
to our data to extract the second order response, and the result is 



d2Re(^o 



-0.62, (Nf = 0) 

0.45, {Nf = 2) d^Re(^o 

0.75, {Nf = 3) dfi^ 
0.34, {Nf = 3,/Xs = 0) 



d/i| 

-0.41, {NF = 3,fis = 0) 



-1.42, {Nf = 0) 
-0.31, {Nf = 2) 
0.00, {Nf = 3) 



(4.43) 



for isoscalar and isovector, respectively. In terms of physical parameters 

d^m _ glCFd^Re{Eo) 



(4.44) 



d/i2 2^3^ d/i2 ' 

which gives second derivatives of order Td^m/d/i^ ~ ±0.02, two orders of magnitude 
smaller than the value measured on lattice. While the masses in both perturbative and 
lattice computations are consistent with the free theory result when all possible error 
sources are taken into account, the difference in the second derivatives with respect to 
isoscalar chemical potential is striking. One has to be careful, however, to define the mass 
the same way in both computations before making any comparison. 



45 



4.5.1 Definition of the screening mass at nonzero density 

The correlator Cz in Eq. (j4.12p is a complicated function of z already at leading order, 
as the exact free theory result in Eq. (14.140 shows. The screening masses appear as poles 
in the momentum space Green's functions, but in practice the correlator is measured in 
configuration space, and only exhibits simple exponential behavior at asymptotically large 
distances, the coefficient of this exponential decay being the screening mass. An effective 
z-dependent screening mass can be defined as 

1 

. -~. (4.45) 

which in free theory approaches the screening mass roughly as 1/z. In lattice simulations 
this complication is often removed by measuring the one-dimensional correlations between 
planar sources. On a finite lattice with periodic boundary conditions the correlator then 
behaves as Cid{z) ~ cosh(— mz), which is used in fitting to extract m from data. 

When chemical potentials are turned on, the situation becomes more complicated. The 
definition of the effective mass in Eq. (14.450 is useless because of the oscillations in the 
correlator which cause Cz to periodically go through zero and to negative values, prevent- 
ing us from taking the 2 — > oo limit. At short distances the fall-off is faster than at ;U = 
because of the cosine term cos(2/x2;) « 1 — 2fj?z'^, but this term does not contribute to the 
decay at longer distances. The asymptotic behavior of the correlator is dominated by the 
coefficient of the exponential fall-off, which corresponds to the real part of the momentum 
space pole location. At large distances this coefficient can be extracted from data, but 
usually one cannot do measurements at arbitrarily large separations because of the finite 
lattice size and the exponentially small value of the correlator. 

When determining the mass from short-distance data the potential oscillatory terms 
have to be included in the fit in order to get reliable estimates. For example, the free 
plane-plane correlator at finite chemical potential behaves as ~ cos{p,z) exp(— Mz) with 
fl = fi[ + fij and M = 2ttT. Simply fitting an exponential of the form Cexp(— mz) 
overestimates the fall-off because of the oscillations, giving 

{m-M){m + Mf = fi^{?,m + M) ^ + M2 ~ if^) ' ^^'^^^ 

The derivatives with respect to are also affected by the unfortunate choice of the fitting 
function. For isoscalar chemical potential p, = 2fj,s, and even for the /x^-independent 
screening mass 2'itT of the free theory Eq. (|4.46|) would give 



T 



d//2 



S 



8T 4 



which incidentally is of the same magnitude as the derivatives measured in [80]. In the 
case of isovector chemical potentials p, = and the oscillations vanish at leading order in 
1/z, so this effect disappears for fj^y- 

As discussed above, the actual lattice simulations are carried out by measuring the 
derivatives of correlators at = instead of the full correlator at /x 7^ 0. The expression 
for the second derivative of C{z) = Aexp{—Mz), given that the first derivative vanishes 
at yu = 0, is [81] 

1 d^C(z) Id^A d^M o d^M 

~ ^-r^ = - ^TT^' (4-48) 



C{z) d//2 Adfi'^ d/x2 d/x2 
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where in the last equahty we have substituted A = cos (2/^2;) and the whole expression is 
written in the limit of infinitely long lattice. The oscillations should then show as quadratic 
z-dependence of the second derivative, but at least the data in [81] does not support the 
existence of this kind of term. 

Given that the second order response does not show signs of oscillation in lattice sim- 
ulations, it is not clear if the cosine terms clearly visible at the leading order and in 
one-loop terms survive at the non-perturbative level. However, understanding the po- 
tentially complicated behavior of the correlator at short distances is necessary to extract 
reliable information from data measured at finite distances. In particular, the oscillatory 
terms suggested by the perturbativc calculations should probably be taken into account 
when fitting the correlator at finite chemical potentials. 
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Chapter 5 

Review and outlook 



Dimensional reduction is an approximate method which is perfectly suited for computing 
static observables in weakly coupled high temperature field theory. It is based on the 
observation that when the temperature is large and the coupling small, there is a clear scale 
hierarchy between the thermal, electric and magnetic scales ttT, gT and g^T^ respectively, 
and makes full use of that hierarchy by applying effective field theory techniques. In 
perturbation theory reorganizing the perturbative expansion by resummation of classes 
of diagrams is necessary in order to regulate the infrared divergences, and this can be 
carried out systematically using dimensional reduction, without worrying about double 
counting diagrams. For non-perturbative lattice simulations the dimensionally reduced 
effective theory is appealing because of the analytic treatment of fermions, lower space- 
time dimension and the integration out of scales ~ ttT, allowing for larger lattice spacing. 
For these reasons, dimensional reduction has been extensively used to compute properties 
of both QCD and clcctroweak theory at high temperatures. However, it is not of use near 
the QCD phase transition, where is large and the weak-coupling hierarchy of scales 
disappears. It should be also noted that dimensional reduction explicitly breaks the Z{?>) 
symmetry of degenerate QCD vacua above Tc, and the fluctuations between these states 
become important close to the phase transition. 

In this thesis, we have studied two applications of dimensional reduction to the standard 
model physics. The pressure of the electroweak theory, which was previously undetermined 
beyond the first terms, was computed up to three-loop level, or 0{g^), in [2,3]. Combined 
with the previously known QCD pressure, this enables us to study the pressure of the full 
standard model above the electroweak phase transition. We have shown that the pressure 
of the symmetric phase goes smoothly through the transition point once the small value 
of the Higgs field mass is taken into account in the effective theory. The pressure of 
the whole standard model to this order lies 10-15% below the ideal gas value, but the 
perturbative expansion docs not converge well because of the large values of QCD and top 
quark Yukawa couplings. As a special case we have studied the weakly coupled SU(2) -|- 
fundamental scalar theory, where the convergence is good as expected. 

The correlation lengths of mcsonic operators at high temperatures have been computed 
in [1] and extended to finite chemical potentials in [4]. We have computed the leading 
order (free) correlators in full QCD, and then formulated a dimensionally reduced effective 
theory which includes the lowest fermionic modes in order to compute the first pertur- 
bative correction of order g'^ to masses. The fermionic theory is seen to correspond to a 
nonrelativistic theory for heavy quarks in 2-1-1 dimensions, and draws inspiration from the 



48 



heavy quark effective theories used to compute the properties of quarkonia. The screening 
mass correction we find is smah and positive, and it depends in a comphcated way on the 
chemical potentials and the number of dynamical quarks. The finite density result is only 
in qualitative agreement with lattice results, and we discuss the potential differences in 
measuring the mass at finite distances. 

Dimensional reduction has proved to be an efficient tool for computing bosonic ob- 
servables at high temperatures both in perturbation theory and combined with non- 
perturbative methods. As a prominent example, the perturbative expansion of the QCD 
pressure has been computed as far as possible by purely perturbative means, and it shows 
good agreement with 4d lattice simulations if the unknown term related to magnetic 
scales turns out to have the right value. This is in accordance with the conjecture that 
the truncated perturbative series is close to the actual result once all dynamical scales 
have entered the computation [82] . Recently, first steps toward computing the missing 
correction have been taken by evaluating the four-loop pressure of massless 0{N) scalar 
theory, requiring for the first time the computation of four- loop sum-integrals [83]. In 
the weak coupling regime dimensional reduction also extends to large values of n/T, for 
the actual requirement is simply that vrT should be larger than other mass scales, which 
for fj. > ttT translates into vrT S> rriE ~ gfJ. as discussed in [84]. An often overlooked 
part of the perturbative computations are finite quark masses [85], which give sizeable 
contributions near the thresholds where they enter. 

New directions for effective theory computations are provided by the possibility of study- 
ing also fermionic operators, as discussed in chapter [H To estimate the reliability of this 
method it would be central to understand the origin of differences between the perturbative 
and lattice results for the masses, possibly also implementing NRQCD3 on lattice. 
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